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FOREWORD 


The  present  study  was  undertaken  with  the  purpose  of  developing  a  local 
continuum  theory  of  fracture  of  sufficient  breadth  to  encompass  the  constitutive 
response  of  degrading  materials  in  the  presence  of  three-dimensional  stress  and  strain 
histories.  Two  specific  materials  were  investigated:  grey  cast  iron  and  plane  concrete 
of  low  to  intermediate  strength. 

Because  continuum  damage  theory  is  a  new  and  emerging  field  and  because  a 
great  deal  of  what  can  be  done  with  such  an  approach  is  still  unknown,  it  was 
thought  wise  to  test  the  fundamental  concepts  (specifically  the  form  of  the  free  energy 
in  the  presence  of  damage)  on  a  brittle  material,  in  our  case  cast  iron,  whose 
constitutive  equation  may  be  considered,  quite  adequately,  elastic.  The  lessons  learned 
from  this  first  study  would  then  be  useful  in  the  more  complicated  task  of  developing 
such  a  theory  for  a  plastic  material  such  as  concrete  where  the  interaction  between 
plasticity  and  damage  was  expected  to  be  much  more  complex. 

This  approach  proved  very  fruitful.  The  work  on  cast  iron  showed  that  the  free 
energy  representation  was  indeed  proper,  in  view  of  the  excellent  agreement  between 
the  calculated  and  experimental  results  in  elastic  domains  with  a  sharp  geometric 
discontinuity  (a  crack).  Since  in  the  endochronic  theory  of  plasticity  the  free  energy 
is  the  basis  for  the  derivation  of  the  constitutive  equation,  this  result  was 
instrumental  in  the  development  of  the  damage  sensitive  constitutive  theory  for 
concrete. 

Within  the  limited  scope  of  the  experiments  performed  the  theory  proved  effective 
in  depicting  the  observed  response  for  both  cast  iron  and  concrete.  Of  course  this  is 
only  a  beginning  and  many  questions  remain.  These  are  discussed  in  somewhat 
greater  detail  in  the  sections  on  accomplished  objectives  and  recommended  further 
work. 


SUMMARY  OF  ACCOMPLISHED  OBJECTIVES 


•  A  continuum  theory  of  damage  first  developed  by  the  author  in  1987  —  in  terms 
of  a  damage  coordinate  and  an  integrity  tensor  —  has  been  refined  and  extended  to 
apply  to  the  constitutive  response  of  brittle  materials  in  the  presence  of  evolving 
damage  consisting  of  sharp  cracks.  The  theory  has  been  tested  by  experiments  on 
grey  cast  iron  thin  plates  with  a  central  crack  normal  and  at  various  angles  to  the 
direction  of  pull.  Calculated  and  measured  fracture  stresses  were  found  to  be  in  close 
agreement  generally,  and  in  reasonably  good  agreement  for  angles  of  inclination  less 
than  30  degrees  to  the  direction  of  pull.  A  finite  element  program  has  been 
developed  for  the  purpose  of  performing  the  calculations. 

•  Of  great  significance  is  the  fact  that  a  conceptual  boundary  has  been  established 
between  the  present  damage  theory  and  Linear  Elastic  Fracture  Mechanics,  in  the 
sense  that  both  theories  give  rise  to  t‘-e  same  result  under  Mode  I  conditions,  when 
the  damage  coordinate  is  the  crack  length. 

•  A  central  key  to  the  success  of  the  theory  is  the  evolution  equation  derived  from 
irreversible  thermodynamics.  Despite  its  continuum  nature  it  predicts  damage  that  is 
very  localized  and  is  confined  to  the  immediate  neighborhood  of  the  crack  tip  in 
agreement  with  observation. 

•  A  micromechanical  connection  has  been  established  between  the  integrity  tensor 
and  the  crack  length  under  Mode  I  conditions.  The  elements  exist  for  the  extension 
of  this  relation  to  more  general  cases. 

•  A  constitutive  theory  of  damage  in  plane  concrete  was  developed  leading  to  two 
different  constitutive  formulations  --  Theories  I  and  II  —  depending  on  the  underlying 
specific  physical  and  geometric  hypotheses.  Theory  II  was  used  to  depict  the  behavior 
of  concrete  in  the  course  of  axial  compression.  Calculated  results  in  terms  of  the 
axial  stress-strain  relation  and  lateral  expansion  were  compared  with  their  experimental 
counterparts.  Close  agreement  was  demonstrated  between  theory  and  experiment. 

•  The  study  yielded  insights  into  the  physics  of  degradation  with  strong  evidence 
that: 


(a)  Substantial  damage  does  not  occur  until  the  material  has  reached 
its  fully  plastic  state. 


(b)  The  rate  of  degradation  is  a  strong  increasing  function  of  the  initial 
elastic  modulus.  This  phenomenon  becomes  particularly  pronounced 
in  the  softening  region  of  the  axial  stress-strain  curve. 


Finding  (b)  has  ramifications  in  regard  to  the  design  of  concrete  structures 
intended  to  survive  damage.  The  inference  is  that  high  strength  concrete  will  degrade 
much  faster  in  the  softening  region  following  the  onset  of  maximum  stress.  This 
point  is  worth  pondering. 

•  An  earlier  damage  evolution  equation  (Valanis,  1987)  proved  very  effective  in 
describing  the  experimental  findings  in  plane  concrete. 


INDICATED  FUTURE  RESEARCH 


•  Directionality  of  Damage. 

We  highly  recommend  the  continuation  of  an  analytical-cum-experimental  research 
program  on  the  constitutive  theory  of  damage  of  plane  concrete,  under  spatially 
heterogeneous  stress  (and  strain_  conditions.  It  is  important  that  we  test  further 
the  predictive  capability  of  the  evolution  equation  in  regard  to  the  directionality  of 
damage  in  spatially  heterogeneous  strain  fields. 


•  The  Constitutive  Representation  of  Damage  in  Situ. 

•  The  question  of  representation  of  observed  damage  by  a  constitutive  damage 

measure  is  still  an  open  question.  In  our  research  we  have  been  able  to  establish 
such  a  relation  in  the  case  of  a  single  crack.  However  how  a  crack  configuration, 
even  in  two  dimensions,  is  to  represented  by  a  damage  measure,  be  it  an  integrity  or 
a  damage  tensor,  is  still  an  open  question.  Much  research  is  needed  in  this  area. 

•  Absorption  of  Transient  Energy  in  Concrete. 

The  finding  that  the  rate  of  degradation  is  a  strong  increasing  function  of  the 
initial  modulus  --  and  by  implication,  the  strength  --  raises  interesting  questions 
I#  regarding  the  absorption  of  energy  in  concrete  due  to  impact.  If  impact  energy  can 

move  away  from  degrading  regions  prior  to  their  disintegration  then  the  absorption 
capability  of  the  medium  will  be  enhanced.  The  interplay  between  the  speed  of 
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propagation  and  the  rate  of  degradation  merits  investigation.  Experimental  (weight 
drop)  and  analytical  studies  on  beams  would  be  a  propos  projects  for  investigating 
this  phenomenon. 

•  The  Ill-posed  Boundary  and  Initial  Value  Problem. 

Spatial  ellipticity  and  temporal  hyperbolicity  of  the  equations  of  the  boundary 
and  initial  value  problem  respectively  is  a  fundamental  issue  in  constitutive  theory.  It 
is  well  known  that  continuum  (local)  damage  theories  without  scale  sensitive  material 
functions  (or  constants)  give  rise  to  ill-posed  boundary  and  initial  value  problems. 
This  question  can  be  addressed  by  the  development  of  a  non-local  damage  theory, 
which  under  homogeneous  conditions  will  give  rise  to  degenerate  forms  with  size 
sensitive  material  parameters  and  functions.  Such  a  theory  would  lend  an  underlying 
unity  to  the  present  research  efforts  in  damage  theories. 

•  Damage  Theory  of  Composite  Materials. 

This  is  an  indicated  extension  of  the  present  damage  theory  of  plane  concrete  to 
deal  with  reinforced  concrete  using  the  theory  of  mixtures.  Certain  inroads  have  been 
made  in  this  area. 

GENERAL  DISCUSSION  ON  THE  QUESTION  OF  DAMAGE 
REPRESENTATION. 

In  recent  years  a  great  deal  of  effort  has  been  devoted  to  the  question  of 
continuum  representation  of  discreet  damage  consisting  of  configurations  of  cracks  and 
voids.  While  this  effort  will  contribute  to  the  understanding  of  damage  there  are 
questions  regarding  the  role  of  such  a  representation  in  a  predictive  theory.  The 
following  discussion  is  an  expression  of  our  views  on  this  matter. 

Let  us  assume  (as  a  matter  of  hypothesis)  that  a  relation  between  a  damage 
pattern  and  a  damage  tensor  has  been  found,  without  using  the  equation  of  evolution 
of  damage.  In  this  event  the  following  situation  would  exist  in  our  case. 


*  Valanis,  K.  C.,  "Thermomechanical  Behavior  of  Anisotropic  Inelastic  Composites. 
A  Micromechanical  Theory,"  J.  Eng.  Mat.  and  Techn.,  To  be  published. 
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Visual  Image  of  Damage  Pattern 


where  $  is  the  integrity  tensor  —  the  damage  measure  used  in  this  work.  Thus  given 
a  visual  image  of  damage  in  a  region,  the  damage  tensor  in  that  region  would  be 
calculated  from  an  analytical  formula. 

Consider  now  the  case  where  the  surface  of  the  domain  has  been  subjected  to  a 
displacement  and/or  traction  history.  Then  at  different  points  along  the  history  the 
domain  can  be  dissected  and  examined  visually  and  can  be  determined  as  a  function 
of  the  path. 

If  the  resulting  strain  field  is  heterogeneous  then  a  sufficient  number  of  sub- 
domains  must  be  dissected  and  visually  examined  to  determine  the  damage  field  in 
the  domain  of  interest. 

To  use  this  information  in  a  predictive  manner,  given  a  strain  history  of  the 
entire  domain  —  we  omit  the  difficulties  of  determining  the  strain  field  in  the  interior 
of  a  domain  —  one  needs  a  constitutive  equation  that  describes  the  material  response 
in  the  presence  of  We  write  such  an  equation  symbolically  as  follows: 


A  "region"  by  definition  will  denote  a  small  sub-domain  of  a  large  domain. 


*  =  F[e;« 


Cm) 

where  F [.]  is  a  function  of  the  history  of  strain  and  0.  Since  both  $  and  e  are 
known  at  every  step  then  a  can  be  calculated  as  a  function  of  the  strain  path. 

If  one  does  not  wish  to  dissect  and  examine  visually  a  large  number  of 
subdomains  of  a  material  domain  and  perform  the  necessary  experiments  to  determine 
the  strain  field  one  must  solve  the  initial  or  boundary  value  problem,  as  the  case  may 
be.  For  this  an  evolution  equation  for  damage  is  necessary.  This  could  be  of  an 
incremental  or  of  functional  type,  but  symbolically  it  may  be  written  in  the  form: 

^  =  G[e;^]  (i  i  i) 

where  G  is  a  function  of  the  history  of  e  and 

It  is  apparent  that  in  the  absence  of  copious  experimentation  and  dissection,  the 
Equivalence  Relation  (i)  is  of  little  value  as  a  predictive  tool  though  it  may  serve  as 
a  proof  test  for  the  evolution  equation  (iii)  by  giving  values  of  ^  (along  the  strain 
path)  which  can  be  compared  with  those  calculated  from  relation  (i).  However  we 
point  out  that  the  constitutive  equation  may  also  be  used  as  a  (more  indirect)  proof 
test  by  comparing  the  calculated  and  measured  stress  response. 

Still  one  can  argue  that  the  relation  (i)  may  be  useful  in  determining  an 
incremental  form  of  the  evolution  equation.  For  instance  one  may  write: 

d f  =  G(f,e)  (iv) 

Knowing  d^  along  a  strain  path  as  well  as  de,  finding  G  would  appear  to  be  a 
plausible  task.  We  explored  such  an  approach  in  the  case  of  cast  iron  but  the 
physics  of  the  problem  was  not  co-operative.  Very  little  observable  damage  was 
present  until  failure.  There  was  no  transient  crack  extension  that  could  be 
measurable  and  thus  be  of  value  in  the  above  approach.  Of  course  our  tests  were 
monotonic.  Cyclic  tests  would  be  more  conducive  to  crack  growth  but  much  less  so 
in  brittle  materials. 

In  ductile  materials  the  situation  is  different.  Cracks  do  propagate  under  cyclic 
strain  histories.  Authors  have  attempted  to  determine  a  relation  between  a  damage 
measure  and  crack  length,  in  Mode  I  conditions,  in  terms  of  the  fracture  toughness 
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parameter  K|c.  Yet  as  is  well  known  and  implicit  in  the  Paris  equation  for  cyclic 
crack  growth  in  metals,  the  rate  of  crack  growth  is  independent  of  K  except  for 
values  of  Kj  nearly  equal  to  K(  .  Thus  the  evolution  of  damage  in  metals,  in  the 
course  of  a  cyclic  history,  bears  no  relation  to  the  damage  measure  determined  from 
monotonic  loading  conditions  involving  . 

Another  instance  where  relation  (i)  would  appear  to  be  of  value  is  the  case  of 
concrete.  Here  damage  does  evolve  in  a  manner  that  is  visible  —  though  this  does 
not  appear  to  be  the  case  under  conditions  of  purely  hydrostatic  pressure  where 
reduction  in  strength  has  been  found  in  the  absence  of  visible  damage.  A  study  by 
Valanis  and  Read  (1989)  and  Valanis  (reported  in  a  proposal  to  AFOSR)  of  the 
experiments  by  Gopalaratnam  and  Shah  (referenced  in  this  report)  and  our  own 
experiments  reported  here  shows  that  the  agent  of  damage  is  not  the  total  strain  but 
the  plastic  strain.  This  must  be  calculated  (in  three  dimensions)  in  the  presence  of 
degrading  elastic  moduli  and  stress  softening.  Thus  the  calculation  calls  for 
knowledge  of  Eqs.  (ii)  and  (iii).  Thus  the  problem  is  coupled  in  the  sense  that  one 
cannot  use  the  equivalence  relation  (i)  to  determine  the  evolution  equation  in  concrete 
without  a  priori  knowledge  of  this  equation. 

It  seems  that  if  one  accepts  the  position  that  the  state  of  damage  is  one  of 
geometry  rather  than  physics,  and  wishes  to  determine  a  continuum  representation  of 
damage,  then  the  approach  indicated  is  to  start  with  a  well-defined  damage  pattern  in 
an  elastic  domain  and  do  tests  to  determine  the  ’’continuum”  average  elastic  response 
of  the  region  in  terms  of  the  integrity  tensor  using  Eq.  (2.1c)  of  Report  I.  Or 
one  may  use  a  finite  element  analysis.  This  is  the  approach  proposed  to  AFOSR. 
One  may  thus  establish  a  relation  between  ^  and  the  parameters  of  a  specific  pattern. 

This  would  be  a  means  cf  resolving  the  important  question  of  initial  damage.  If 
a  domain  is  initially  damaged,  the  value  of  must  be  known  before  the  initial  or 
boundary  value  problem  can  be  solved. 

The  Question  of  Element  Size. 

This  question  has  been  discussed  in  the  literature  a  great  deal  in  the  last 
decade.  What  is  often  left  unsaid,  or  not  realized,  is  that  there  are  two  questions  of 


size.  One  involves  the  effect  of  the  size  of  the  domain  on  the  representation  of 
damage.  Obviously  such  an  effect  exists  and  the  only  question  is  how  to  calculate  it. 
The  effect  will  depend  on  the  damage  pattern.  Thus  if  damage  is  uniformly 
distributed,  the  effect  will  be  small.  On  the  other  hand  if  it  is  concentrated,  the 
effect  will  be  large.  The  damage  pattern  is  a  material  property  and  will  vary  from 
one  material  to  the  other  given  the  same  strain  history.  Thus  the  size  effect  is  a 
geometric  as  well  as  a  constitutive  property. 

The  other  size  effect  concerns  the  computation  of  the  material  response  by 
means  of  a  finite  element  method  in  the  presence  of  damage-induced  softening.  The 
numerical  solution  of  the  associated  initial  or  boundary  value  problem  is  sensitive  to 
the  size  of  the  mesh.  Local  continuum  damage  theories  are  susceptible  to  this  effect 
and  give  rise  to  ill-posed  boundary  or  initial  value  problems. 

Thus  accounting  for  the  effect  of  the  size  of  the  domain  in  the  continuum 
representation  of  damage  does  not  eliminate  the  effect  of  mesh  size  in  a  finite  element 
analysis  if  the  damage  theory  is  of  "local”  character.  In  other  words  local  damage 
theories  apply  only  to  situations  where  the  overall  stress  and  strain  fields  are  spatially 
uniform.  However,  local  theories  are  very  useful  in  understanding  material  behavior. 
They  are  in  fact  special  (degenerate)  cases  of  the  constitutive  response,  when  the 
stress  and  stain  fields  are  spatial  uniform. 

In  conclusion  it  is  our  position  that  what  is  called  for  is  a  non-local  evolution 
equation.  Such  an  equation  would  have  the  effect  of  eliminating  difficulties  associated 
with  the  size  of  the  mesh  in  finite  element  calculations  and  would  resolve  the 
question  of  the  ill-posedness  of  the  initial  and  boundary  value  problems.  It  is  our 
suspicion  that  such  an  equation  will  also  resolve,  in  passing,  the  question  of  the  size 
of  the  domain  on  the  continuum  representation  of  damage. 

We  are  working  on  such  an  approach  right  now  and  plan  to  make  this  the  focus 
of  our  future  work  on  damage  theory.  Our  initial  results  are  very  encouraging. 


A  THEORY  OF  DAMAGE  IN  BRITTLE  MATERIALS 


ABSTRACT 


The  present  paper  addresses  the  geometric  representation  of  damage  in  brittle 
solids,  its  equation  of  evolution  and  its  incorporation  in  the  constitutive  equation. 
The  concept  of  a  damage  coordinate  is  introduced  and  a  thermodynamic  derivation  of 
the  evolution  equation  follows.  A  continuum  damage  theory  ensues.  The  theory  is 
applied  to  the  case  of  a  thin  plate  with  a  central  crack.  Two  cases  are  treated. 

In  Case  (a)  the  tensile  load  is  in  a  direction  normal  to  the  plane  of  the  crack. 
The  fracture  stress  is  calculated  numerically  using  a  finite  element  code  for  various 
crack  lengths  and  a  comparison  is  made  with  its  observed  values.  Excellent 
agreement  between  the  two  is  obtained  with  the  aid  of  only  one  damage  parameter, 
the  damage  propensity  constant.  A  congruence  between  the  theory  and  linear  fracture 
mechanics  in  Mode  I  is  then  shown  when  the  damage  coordinate  is  the  crack  length. 

In  Case  (b)  the  crack  length  is  constant  but  the  angle  of  inclination  of  the  crack 
to  the  direction  of  pull  is  varied.  The  fracture  stress  is  calculated  as  above  for 
various  angles  and  comparison  is  made  with  observed  values.  While  overall  agreement 
is  reasonable,  excellent  agreement  is  found  when  the  angle  is  greater  than  45*. 
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1.  INTRODUCTION 


Damage  plays  an  essential  role  in  the  constitutive  behavior  of  materials  Its 
nature,  in  terms  of  its  measurement  and  mathematical  representation  is  complex.  It 
can  vary  from  a  well  defined  crack  to  randomly  distributed  damage  that  can  be 
thought  of  as  a  ’’damage  field”.  It  is  important  to  understand  under  what  initial 
damage  conditions,  stress  history  and  material  properties,  one  type  will  prevail  over 
the  other,  or  when  the  two  types  will  co-exist. 

Constitutive  theories  in  the  context  of  continuum  damage  have  been  proposed  by 
Dragon  and  Mroz  [1],  Lemaitre  and  Mazars  [2],  Chaboche  [3],  Krajcinovic  and 
Fouseka  [4]  and  Fouseka  and  Krajcinovic  [5]  and  Krajcinovic  [6].  In  this  regard  see 
also  a  recent  paper  by  Chaboche  [7]. 

Since,  however,  failure  often  occurs  by  the  extension  of  one  dominant  crack,  at 
times  in  the  presence  of  distributed  damage,  a  comprehensive  theory  must  be  able  to 
deal  with  both  types,  and  their  mixture,  simultaneously.  Specifically  a  continuum 
theory  should  be  able  to  deal  with  damage  due  to  well-defined  sharp  cracks.  The 
theory  in  this  study  is  being  pursued  with  this  goal  in  mind.  The  state  of  directional 
damage  (in  the  form  of  oriented  cracks)  is  represented  by  an  integrity  tensor  g,  and 
random ly  distributed  scalar  damage  by  a  scalar  function  fl,  both  of  which  may  co¬ 
exist. 


While  both  types  of  damage  were  dealt  with  in  a  previous  work  Valanis  [8],  in 
this  study,  which  is  reported  in  Part  I  of  this  Report,  we  concentrate  on  damage  in 
brittle  solids,  which  is  directed,  and  begin  with  the  case  of  a  well-defined  crack.  The 
reason  for  this  is  born  out  of  the  foregoing  discussion.  In  addition,  however,  there  is 
the  question  of  Linear  Fracture  Mechanics  (LEFM)  which  has  been  successful  in 
dealing  with  fracture  due  to  sharp  cracks  in  brittle  materials  and  in  Mode  I 
conditions.  In  developing  the  theory  we  feel  that  we  must  establish  a  conceptual 
boundary  with  LEFM. 

To  this  end  we  use  the  continuum  damage  theory  developed  herein,  to  solve  the 
problem  of  a  plate  under  tension  with  a  central  crack.  Two  case  are  considered. 

Case  (a). 

In  this  case  the  crack  is  normal  to  the  direction  of  pull.  The  fracture  stress  is 
determined  as  a  function  of  the  crack  length.  Its  values  are  then  compared  with 
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those  obtained  by  experiment  on  grey  cast  iron,  and  from  the  analytical  linear  elastic 
fracture  mechanics  solution.  The  agreement  between  the  damage  theory  and 
experiment  is  excellent.  The  fracture  mechanics  solution  disagrees  with  the  experiment 
for  values  of  the  crack  length  that  are  small  but  is  otherwise  in  close  agreement  also. 

Case  (b). 

In  this  case  the  crack  is  of  fixed  length  but  its  angle  of  inclination  to  the 
direction  of  pull  is  varied.  The  agreement  between  the  damage  theory  and  experiment 
is  good  for  angles  of  inclination  that  exceed  45  degrees.  A  numerical  solution  could 
not  be  obtained  for  angles  smaller  than  30  degrees  because  of  slow  convergence  of 
the  iterative  algorithm.  Since  the  solution  for  an  angle  of  zero  degrees  is  known,  this 
being  the  solution  for  the  plate  without  a  crack,  an  approximate  analytic  solution  for 
angles  less  than  30  degrees  has  been  obtained  by  interpolation.  A  deviation  between 
theory  and  experiment,  however,  is  already  apparent  for  values  of  inclination  between 
45  and  30  degrees.  While  the  reason  for  this  is  discussed  at  greater  length  in 
Section  7  we  feel  that  the  primary  cause  is  the  finiteness  of  the  crack  width  and  the 
non-zero  value  of  the  radius  of  the  crack  tip.  This  could  not  be  avoided  as  the  crack 
was  made  with  best  available  technology.  See  Appendix  B. 

In  addition  to  these  overall  results,  the  theoretical  study  has  brought  forth  two 
central  findings  which  are  the  following: 

(i)  The  evolution  equation  —  with  a  firm  basis  in  thermodynamics  in 
the  sense  that  it  abides  by  the  law  that  the  rate  of  change  of  a 
thermodynamic  variable  is  proportional  to  the  dual  thermodynamic 
force  —  gives,  in  the  course  of  increasing  tension,  a  damage  field 
which  is  very  localized  and  is  confined  essentially  to  the  vicinity  of 
the  crack  root,  in  accord  with  observations  on  brittle  fracture  and 
our  own  experiments  on  grey  cast  iron.  Visual  examination  of  the 
material  with  a  microscope  --  magnification  400  --  in  the 
neighborhood  of  the  crack  root,  after  fracture,  revealed  no  visible 
changes  in  material  structure,  except  very  close  --  within  a 
millimeter  or  so  --  of  the  crack  root. 
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We  have  therefore  a  conceptual  as  well  as  a  quantitative  confluence  among 
(a)  the  continuum  theory,  (b)  Linear  Elastic  Fracture  Mechanics  and  (c)  experimental 
observation  in  Mode  I  conditions.  We  also  have  good  agreement  between  the  damage 
theory  and  experiment  in  the  case  of  inclined  cracks  for  angles  of  inclination  greater 
than  45  degrees.  The  conclusion  on  smaller  angles  must  await  further  developments. 


2.  DAMAGE  THEORY. 


We  begin  by  setting  out  the  fundamentals  of  a  theory  for  incorporating  directly 
the  effects  of  damage  in  the  constitutive  equation  of  a  linear  elastic  material.  A  work 
which  addresses  the  behavior  of  elastic-fracturing  and  plastic-fracturing  solids  has 
appeared  in  two  unpublished  reports  by  Valanis  [8j  and  Valanis  and  Read  [9J.  Here 
we  shall  be  concerned  only  with  directed  damage.  The  question  of  scalar  damage  in 
the  presence  of  directed  damage  was  dealt  with  previously  by  Valanis  [8]. 

The  fundamental  premise  of  the  theory  is  that  directed  damage  may  be 
represented  by  a  tensor  valued  function  "the  integrity  tensor",  which  appears 
explicitly  in  the  equation  for  the  free  energy  and  the  constitutive  equation. 
Thermodynamically,  g  is  an  internal  variable.  In  the  case  of  elastic-fracturing  solids 
the  free  energy  density  is  given  by  Eq.  (2.1a)  and  the  constitutive  relation  by  Eq. 
(2.1c).  Specifically  [8], 

r  ■  WQtijtMfij'u.  *  ikVijV  <21a> 


where  X  and  n  are  the  Lame  constants.  Thus,  since 


it  fol lows  that: 

ffij  =  X*ijVek£  +  2^ik*j  l€kt 


(2.1b) 


(2.1c) 


Equation  (2.1a)  is  a  statement  to  the  fact  that 
*  =  2  Cijk£  eij  ek£. 


(2.2) 


where  C,  the  stiffness  tensor,  is  homogeneous  and  quadratic  in  g. 

There  are  two  fundamental  questions  regarding  the  tensor 

(i)  What  is  its  geometric  nature?  More  specifically,  how  is  it  related 
to  the  actual  material  damage  at  some  microscale  —  i.e.,  a  scale 
smaller  than  the  one  that  underlies  the  macroscopic  constitutive 
behavior  of  the  material? 


(ii)  How  does  $  evolve  with  the  history  of  deformation? 


In  regard  to  question  (i).  it  was  previously  shown  Valanis  [8],  that  if  in  a 
material  element  a  principal  value  of  g,  say  $  ,  is  zero,  and  INI  is  the  eigenvector  of 
<t>  .  then  the  traction  on  a  surface  normal  to  is  zero.  That  fs.  the  element  cannot 

'  .  r 

support  shear  or  direct  stress  on  that  surface.  Thus,  $  represents  damage  on  a 
plane  normal  to  j\L  and  as  such  it  is  a  measure  of  the  effective  area  in. the  presence 
of  microcrack(s)  in  that  plane.  Thus  when  ^  =  0,  an  effective  plane  crack  has 
developed  across  the  entire  element  on  a  plane  normal  to  (VI  .  Conversely  if  <p  =  1 
no  damage  has  occurred  on  that  plane. 

With  regard  to  question  (ii),  i.e.,  the  evolution  of  g,  the  following  equation  was 
proposed  as  a  basis  for  the  description  of  damage  due  to  mechanisms  of  brittle  (or 
semi-brittle)  fracture: 

=  *  [q“)'V  (2-3) 


The  notation  in  Eq.  (2.3)  is  the  following:  dfSa  is  the  change  in  principal  component 
of  g  due  to  a  strain  increment  dg  and  n?  are  the  eigenvectors  of  d£.  Also,  is 
a  damage  coordinate  and  Q®  a  "damage  force"  in  the  direction  na,  i.e.. 


na  a  an 

Q  =  n.n.Q. . 
"  i  j^ij 


(r  not  summed) 


(2.4) 


The  power  index  m  in  Eq.  (2.3)  is  a  material  constant. 
The  damage  coordinate  {a  is  given  by  Eq  (2.5)  i.e., 


fkdea  ;  dea>0,  ea>0,  Qa>0 
'  n  '  n  '  m 

.0  ,  otherwise 


(2.5) 


where  =  e.-n^n*?  (a,  not  summed)  and  k  is  a  positive  scalar.  A  complete 

discussion  of  these  equations  when  Q  =  g  is  given  by  Valanis  (8j. 


Further  Discussion  of  Eq.  (2.2). 

The  appropriate  "driving  force"  is  of  central  importance  to  the  evolution 
equation.  In  Ref.  8,  Q  was  set  equal  to  g  because  the  resulting  relation  was  of  a 
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type  that  governs  annihilation  of  species  and,  therefore,  degradation.  In  brittle 
materials  such  as  grey  case  iron,  however,  and  in  the  case  where  damage  is  due  to  a 
well-defined  crack,  we  found  that  this  equation  overestimates  the  damage  in  zones  of 
low  strain  intensity  [8].  This  result  is  contrary  to  the  behavior  of  brittle  materials 
and  our  own  observations  on  grey  cast  iron. 

The  Physical  Fracturing  Process. 

The  process  of  fracture  of  brittle  and  semi-brittle  materials,  such  as  grey  cast 
iron  which  exhibits  a  small  amount  of  plasticity  but  fractures  at  very  low  strain 
(~  0.5%),  shows  two  essential  characteristics  which  are  always  consistently  and  clearly 
observable: 

(i)  The  damage  is  at  first  very  minimal,  but  then  increases  suddenly 
becoming  catastrophic. 

(ii)  In  the  case  of  a  single  crack,  material  regions  away  from  the  crack 
remain  essentially  damage-free  during  the  deformation  process  while 
damage,  which  leads  to  catastrophic  failure,  remains  confined  to  the 
tip  of  the  crack. 

The  above  characteristics  must  be  depicted  quantitatively  by  the  damage 
evolution  equation  whose  suitability  or  otherwise,  will  be  determined,  in  part,  by  its 
capability  to  describe  the  essential  physics  of  the  damage  process. 

A  Thermodynamic  Approach. 

To  arrive  at  a  suitable  evolution  equation  we  appeal  directly  to  thermodynamics. 
We  recall  that  since  £  is  an  internal  thermodynamic  variable,  the  thermodynamic  force 
Q  that  drives  the  damage  process  is  given  by  the  relation 

9  ■  -  tjj  (2-6) 


Also  the  direction  of  the  damage  increment  -  dg  is  dictated  by  the  eigendirections  na 
of  the  increment  of  strain  d£  In  view  of  the  above  observations  we  make  the 
following  two  fundamental  stipulations: 

(i)  dg  and  dg  are  coaxial. 
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(ii)  The  operative  force  drivmg  the  damage  process  in  the  direction  nQ  is  Q^. 
which  is  Q  projected  in  direction  na  in  accordance  with  Eq.  (2.3).  Using  then  the 
basic  linear  thermodynamic  law  "that  the  rate  of  change  of  a  thermodynamic  variable 
is  proportional  to  the  conjugate  thermodynamic  force,”  we  obtain  the  relation: 


d^ 


(2.7) 


where  b  is  a  material  scalar  and 


(2.8) 


k  having  the  significance  of  a  damage  propensity  parameter.  As  before.  dea  are  the 
eigenvalues  of  de. 

Previously,  when  Eq.  (2.2)  was  used,  Valanis  [8],  the  following  were  the 
conditions  that  d£a  be  non  zero: 


de“  >  0 


ea  = 
n 


e. -nana  >  0 
•J  '  J 


(2.9) 


Otherwise,  d£a  =  0.  However  in  the  present  case  the  explicit  appearance  of  the 
thermodynamic  force  in  the  evolution  equation  brings  the  Clausius-Duhem  dissipation 
inequality,  i.e.. 


di  >  o 


\m  * 0 


(2.10) 


into  play  directly.  Thus  in  view  of  Eqs.  (2.6)  and  (2.7)  and  the  fact  that 


d^. .  =  y  d^an?n? 
r'J  ~  r  i  J 


(2.11) 


•♦a 

n 


being  the  eigenvectors  of  dg.  it  follows  that 
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d*  =  - 

a 


a  a  .  A 
g  —  n .  n .  y  0 
09. ■  I  J 

f'J 


(2.12) 
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Since  ineq.  (2.12)  must  be  true  for  each  a  individually,  i.e.,  for  each  independent 
fracture  mechanism,  then  the  following  inequality  must  be  satisfied  for  all  a: 

*  ntnt  >  0  (2-13) 

Yi  j  J 

Furthermore,  since  d$a  is  always  negative  whenever  it  is  nonzero,  then 

4fr  > 0  <2 14> 

for  all  a,  whenever  |d^a|  t  0. 

We  have  thus  obtained  the  full  set  of  conditions  that  govern  the  evolution  of 
the  damage  coordinate  d£a.  Thus: 

df“  =  kde°  (2.15) 

whenever : 

dea  >  0  ,  e*  >  0  ,  n*n®  >  0  (2.16) 

rij  J 


but  dfa  =  0  otherwise. 

We  recall  that  the  search  for  a  new  evolution  equation  was  motivated  by  the 
fact  Eq.  (2.2)  overestimated  the  rate  of  growth  of  damage  in  regions  of  low  strain 
intensity.  To  see  that  Eq.  (2.7)  is  superior  to  Eq.  (2.3)  in  this  respect  we  write  it 
in  its  explicit  form  given  below: 

V  =  -  bdf“{x£ij  (^k£J  .  2M^ikej}"°n“  (2.17) 

Note  that  the  righthand  side  of  Eq.  (2.17)  is  now  of  order  ||e||  --  whereas  previously 
it  was  of  order  |ig||°.  Other  things  being  equal  the  rate  of  diminution  of  $a  will, 
therefore,  be  expected  to  be  substantially  higher  in  regions  of  high  strain  intensity 
such  as  the  root  of  an  existing  crack. 

Material  Parameters  and  Their  Experimental  Determination. 

It  may  be  seen  from  Eq.  (2.7)  and  (2.8)  that  b  and  k  merge  into  a  constant  bk 
so  that  without  loss  of  generality  we  take  b  to  be  equal  to  unity.  There  are, 
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therefore,  three  constants  to  be  determined:  k,  the  fracture  susceptibility,  and  the 
two  elastic  constants  X  and  /*.  In  a  uniaxial  test  Young’s  modulus  E  was 
determined,  while  u,  the  Poisson  ratio,  was  taken  equal  to  0.3.  The  elastic  constants 
X  and  fi  were  then  calculated  and  were  found  to  be  15  x  10  and  10  x  10  ksi, 
respectively. 

The  constant  k  is  determined  from  a  simple  tension  test  on  an  uncracked 
specimen.  Under  these  conditions 


a^  =  Eei* 


lrl 


(2.18) 


and  using  Eqs.  (2.17)  and  (2.18) 


=  kEe^dej 


(2.19) 


Equation  (2.19)  ma  be  integrated  to  give: 


*  is  e 


-(kE/3)e: 


(2.20) 


I 


I 


Thus,  in  view  of  Eqs.  (2.18)  and  (2.20) 
al  =  Ee1  exp[(2k/3)EeJ] 

Fracture  will  take  place  when  a  reaches  a  maximum  value. 


(2.21) 

The  condition  is 

(2.22) 


which,  in  conjunction  with  Eq.  (2.21)  gives  the  relation: 

£1  «„■  d/2«E>1/3  <2'23> 

"l  *  Ed/2k.E)1/3  (2-2") 

where  e  in  Eq.  (2.24)  is  the  base  of  the  natural  logarithm. 
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Two  different  batches  of  specimens  were  used  in  the  experiments,  one  from 
ingot  (a)  for  Case  (a)  and  one  from  ingot  (b)  for  Case  (b). 

Case  (a>. 

In  a  series  of  experiments  grey  cast  iron  specimens  from  ingot  (a)  were  tested 

in  axial  tension  and  the  fracture  stress  a,  was  found  to  lie  within  the  limits 

1  max 

32  x  103  <  a .  <  38  x  103  psi 

1  max 

while  an  upper  bound  for  E  was  found  to  be  26.6  x  106  psi.  Thus,  choosing  k  to 
be  equal  to  2000  in  units  of  (ksi)’1,  a  was  found  to  be  40.2  ksi  which  is  close 
to  the  upper  bound  of  the  experimentally  determined  value. 

Case  (b). 

From  tensile  experiments  on  specimens  from  ingot  (b).  E  was  found  to  be 

17  x  106  psi  while  values  of  the  fracture  stress  clustered  close  to  a  value  of 
3  3.-1. 

30  x  10  psi.  Thus  k  was  found  to  be  equal  to  2  x  10  (ksi)  in  this  case  as  well. 

Values  of  X  and  /<  were  then  calculated  for  both  cases  by  setting  v  =  1/3. 

COMPUTATIONAL  STUDY 

As  discussed  in  the  Introduction,  the  computational  study  consisted  of: 

(a)  Developing  a  Finite  Element  Analysis  Program,  in  the  presence  o*  the 
damage  evolution  equation,  for  the  computation  of  the  deformation,  stress 
and  damage  fields  in  a  flat  plate  in  the  presence  of  a  central  crack. 

(b)  Computing  these  quantities  for  central  cracks  of  various  lengths  in 
Case  (a)  and  for  a  crack  of  constant  length  and  varying  angle  of 
inclination  to  the  direction  of  pull  in  Case.  (b). 


(c)  Developing  a  numerical  Algorithm  for  the  above  purpose. 

Case  (a). 

The  numerical  experiment  performed  was  one  of  the  longitudinal  displacement 
control  (in  the  presence  of  free  transverse  displacement)  on  the  outer  boundary.  All 
other  boundaries  were  stress-free.  An  incremental  displacement  of  2.5  x  10°  units 


was  imposed  at  each  stage.  The  computations  were  carried  out  for  a  square  grid  of 
fixed  size.  For  reasons  relating  to  the  geometry  of  the  specimen.  3/64"  was 
designated  as  a  unit  of  length.  Half-domains  were  used  in  the  computations  utilizing 
the  longitudinal  axis  of  symmetry. 


To  test  the  effect  of  specimen  length  on  the  fracture  stress  two  such  domains 
were  used.  The  first  had  a  width  of  eight  units  and  a  height  of  eleven  units  of 
length.  The  second  had  the  same  width  as  the  first,  but  a  height  of  twenty-one 
units.  Stress,  strain  and  damage  fields  were  obtained  in  both  cases.  For  the 
domains  tested  the  length  effect  was  negligible.  Fracture  was  considered  to  have 
occurred  when  the  displacement  field  was  such  that  two  halves  of  the  plate 
(symmetrically  located  relative  to  the  crack  line)  began  to  move  relative  to  each  other 
as  rigid  bodies. 


Case  (b). 


In  this  case  the  boundary  conditions  were  the  same  as  in  Case  (a).  When  the 
crack  was  inclined,  however,  there  was  no  reflective  axis  of  symmetry.  For  this 
reason  the  entire  plate  (width  16  units,  length  32  units)  was  used  in  the  calculations. 
In  Figure  9  we  show  the  geometry  of  the  right  half  of  the  plate  and  the  "idealized” 
set  of  cracks  used  in  the  computation  --  in  the  sense  that  some  of  the  cracks  were 
slightly  longer  than  the  ones  used  in  the  experiments.  The  resulting  advantage  was 
that  in  the  idealized  configuration  all  elements  affected  by  the  presence  of  a  crack  has 
the  same  initial  value  of  i.e.. 


fc2 


in  ~ 


cs  s 
0  0 


cs  0 
2  „ 


(2.24a) 


where  c  =  cosfl,  s  =  sin0  and  9  =  t/2  -  a,  where  a  is  the  inclination  of  the  crack 
to  the  direction  of  pull.  In  this  case  failure  was  considered  to  have  occurred  when 
the  axial  stress  in  the  element  adjacent  to  the  crack  reached  a  maximum  value,  since 
following  this  event  the  displacement  field  became  prohibitively  large  and  erratic  and 
while  convergence  of  the  iterative  process  failed. 


Brief  Description  of  the  Numerical  Algorithm. 

The  boundary  value  problem  for  a  brittle  material  in  the  presence  of  initiation 
and  evolution  of  damage  is  non-linear,  and  there  are  strong  coupling  effects  between 
the  elastic  stiffness  and  the  state  of  damage.  T.ie  coupling  is  manifested  through  the 
constitutive  matrix  [k  ,]  which  is  directly  dependent  on  the  integrity  tensor  fc. 
Specifically  the  constitutive  response  in  plane  stress  is  given  by  Eq.  (2.25)  where 


°iil  en 

°22  ~  £22 

a12J  l712. 


(2.25) 


where  7  is  the  engineering  shear  strain  and 

4>h&  *  2m)  (x#h#22  *  2^i2)  ^ii^22(X  +  W 

[^]  =  (x^2  -  2 ^22(X  +  2^  ^12(X  +  2^  (2  26) 

^11^12 (X  +  2/i)  ^12^22(X  +  2/i)  (X  ♦  m)^2  + 

On  the  other  hand  the  evolution  equation  for  the  integrity  tensor  £  involves  both  the 
state  of  damage  (given  by  £)  the  state  of  strain  e  as  well  as  the  increment  of  strain 
de  which  determines  the  damage  coordinate,  increments  d£a  (a  =  1,2.3).  Specifically: 

d^Q  =  Q^d£fl  (a  not  summed)  (2.27) 

<  <2-28) 

«ij  *  •  S)  *  2^ik£j£.  (2  «) 


The  algorithm  adopted  here  was  one  of  incremental  application  of  boundary 
displacements  in  conjunction  with  the  Newton-Raphson  method  of  iteration.  In  the 
equilibrium  configuration  each  node  is  in  equilibrium.  Loads  applied  externally  plus 
loads  applied  by  the  elements  must  sum  to  zero.  If  this  sum  if  [AR J  instead  of  the 
null  vector,  the  load  imbalance  [ AR)  will  produce  a  displacement  difference  vector  [ Au] 
which  is  computed  by  solving  the  equation 
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[K0][Au]  =  [AR] 


(2.30) 


where  [KQJ  is  the  stiffness  matrix  of  the  structure  in  absence  of  damage.  In  iterative 


notation: 


*  [K^'hsRlS 


where  k  denotes  the  iteration  number,  n  the  incremental  loading  step  and 


(2.31) 


(2 . 32) 


i.e.,  [ARP  is  calculated  using  the  current  secant  stiffness  matrix  [K  ,]  and  the  current 
iterated  displacement  uk.  ^ 

The  imbalance  "node  force”  AR  and  the  displacement  difference  between  two 
iterations  will  decrease  during  the  iteration  process.  Computation  is  stopped  when  a 
tolerance  is  reached,  i.e.,  when 


i  k+1  k.  ,  .  k. 

u  -  u  <  e  u 
1  n  n  1  1  n 1 


(2 . 33) 


where  e  is  a  small  number.  In  this  study,  e  was  set  equal  to  10  . 


As  discussed  previously,  in  Case  (a)  the  plate  was  considered  as  having 
fractured  when  the  difference  in  the  displacements  of  any  two  nodes  in  either  the  top 
or  bottom  part  of  the  plate  relative  to  crack  line,  was  less  than  10  6.  while  in  Case 
(b)  fracture  was  considered  to  have  occurred  when  the  stress  in  the  element  adjacent 
to  the  crack  reached  a  maximum  value.  The  numerical  computations  were  conducted 
on  a  VAX  750. 

Discussion  of  Results. 

Figure  1  shows  the  calculated  curve  of  fracture  stress  versus  crack  length,  based 
on  the  damage  theory.  Also  shown,  are  the  experimentally  determined  values  of 
fracture  stress.  The  comparison  shows  excellent  agreement  between  calculated  and 
experimental  data.  Figure  la  shows  the  calculated  fracture  curve  when  linear  elastic 
fracture  mechanics  is  used.  It  is  apparent  that  LEFM  gives  good  agreement  with  the 
data  for  larger  values  of  crack  length,  as  expected. 
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Figure  la.  Comparison  of  LEFM  prediction  with  experimental  observation. 
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Figures  2,  3.  and  4a  show,  respectively,  the  stress,  strain  and  damage 
distribution  along  the  crack  line  when  the  half-crack  is  one  unit  long.  Also  shown  in 
Figure  4a  is  the  damage  distribution  along  the  elements  just  below  the  crack  line  — 
again  when  the  half-crack  length  is  one  unit.  It  may  be  seen  that  the  calculated 
damage  away  from  the  crack  is  indeed  negligible.  In  Figure  4b  we  show  the  damage 
distribution  along  the  crack  line  when  the  half-crack  is  one  unit  long  and  three  units 
long  —  for  comparison. 

Figure  5  shows  the  stress-strain  curve  for  the  element  adjacent  to  the  tip  of  the 
crack  when  its  half  length  is  three  units.  Material  softening  is  clearly  visible  leading 
to  failure  when  the  stress  reaches  a  maximum  value.  Figure  6  shows  a  similar  type 
of  behavior,  when  the  stress  of  the  element  adjacent  to  the  crack  tip  is  plotted  versus 

boundary  displacement  --  for  a  crack  whose  half-length  is  again  3  units. 

The  attenuating  effect  of  damage  on  the  stress  concentration  is  shown  clearly  in 
Figure  7.  while  Figure  8  shows  the  evolution  of  damage  in  the  element  adjacent  to 

the  crack  tip  versus  boundary  displacement.  It  may  be  seen  that  damage  proceeds  at 

a  slow  rate  in  the  range  of  small  strains,  accelerating  rapidly  to  a  catastrophic  rate  as 
the  strain  increases.  The  half  length  of  the  crack  was  again  three  units. 

In  Figure  9  we  show  a  comparison  between  the  experimental  and  calculated 
values  of  the  fracture  stress  for  various  angles  of  inclination  a  of  the  crack  to  the 
direction  of  pull.  Agreement  between  theory  and  experiment  is  very  good  for  values 
of  a  greater  than  45  degrees.  We  believe  that  the  discrepancy  for  values  of  a  smaller 
than  45  degrees  is  due  to  the  finiteness  of  the  crack  width  and  root  radius.  This 
conjecture  was  arrived  at  by  observing  that  the  crack  extension  did  not  originate  at 
the  root  of  the  crack  but  at  a  point  where  a  vertical  line  made  tangential  contact 
with  the  crack  contour.  This  question  needs  further  investigation. 


Figure  2.  Stress  distribution  along  crack  line. 


Distance  to  Crack  Tlpfriunit  length} 

Figure  3.  Strain  distribution  along  crack  fine. 


Figure  4.  Damage  distribution  along  and  below  crack  line. 
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Dlstonce  to  Crack  Tip(in  unit  length) 

Figure  4b.  Damage  distribution  along  crack  line. 


Stress-strain  curve  for  the  tip  of  crack. 


Boundary  Displacement  X10  (in  unit  length) 

Figure  6.  Stress  versus  boundary  displacement. 


Figure  7.  Stress  concentration  factor  versus  boundary  displacement. 
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Figure  8.  Damage  evolution  versus  boundary  displacement. 
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Figure  9.  Fracture  stress  versus  angle  of  crack. 


3.  DAMAGE  THEORY  IN  THE  CONTEXT  OF  FRACTURE  MECHANICS  IN 
MODE  I. 


As  pointed  out  in  the  Introduction  our  purpose  in  this  study  is  two-fold.  We 
wish  to  examine  the  proposed  theory  in  the  light  of  experimental  data,  as  well  as  in 
the  context  of  linear  fracture  mechanics  and  thereby  give  perspective  to  the  equations 
of  the  theory.  Moreover  we  are  also  aiming  at  determining  how  £  depends 
quantitatively  on  the  microdamage  when  this  consists  of  a  well-defined  crack  of  a 
certain  extent  and  direction  relative  to  a  material  microelement.  Of  course  a  sub-part 
of  our  micromechanical  aim  is  to  determine  the  effective  crack  "extent”  when  the 
damage  is  not  mathematically  well  defined.  However,  this  last  part  must  await  the 
findings  of  the  first  phase  of  this  work. 

We  begin  by  considering  a  thin  sheet  of  "large"  extent  (and  unit  thickness) 

under  axial  stress  in  the  presence  of  a  through  crack  normal  to  the  direction  of 

stress.  Also  within  the  sheet,  is  a  circular  region  whose  center  is  at  the  midpoint  of 

the  crack  and  whose  radius  is  R.  The  strain  energy  of  the  circular  region  in  the 

absence  of  the  crack  is  $  where 

o 


(3.1) 


According  to  the  Griffith  calculation  [10],  in  the  limit  of  very  large  R,  or  vanishingly 
small  "a”,  the  reduction  A¥  in  strain  energy  due  to  the  onset  of  the  crack  at  fixed 
strain  e  is  given  by  Eq.  (3.2). 


A*  =  ±  jr 


(3.2) 


I 

Hence,  the  free  energy  of  the  circular  region  in  the  presence  of  the  crack  is: 


-  M  .  1  £  »[R2  -  a2] 


(3.3) 


or 

*  “  2  1“  II  "  (a/R)2^2 


(3.4) 
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or 


*  = 


(3-5) 


since  the  far  field  axial  stress  is  related  to  the  far  field  axial  strain  by  the  relation: 
a  —  E  e. 


In  the  Griffith  fracture  theory  the  "energy  release  rate"  at  fracture  is  equal  to 
the  rate  of  creation  of  surface  energy,  i.e.. 


(3.6) 


since  for  an  increment  of  crack  length  "da"  the  surface  energy  created  is  27  da. 
Thus  using  Eqs.  (3.5)  and  (3.6)  in  the  limit  of  R  “, 

4^  =  27  (3.7) 

which  is  the  Griffith  equation  for  the  critical  stress,  i.e.. 


a 

c 


2Ey 
<  ira 


(3.8) 


The  Present  Theory. 


It  is  easily  shown  that  in  the  case  of  a  thin  plate  under  axial  stress  conditions 
=  0.  i  i  j.  $n  =  *33  =  1)  Eq.  (2.1a)  gives: 

*=±Ee2^  (3.9) 


Thus  $  which  is  the  energy  of  a  region  with  radius  R  becomes 


$  =  *R2  f  =.  |  Ee2^2wR2 


(3.10) 


comparison  of  Eq.  (3.10)  with  Eq.  (3.4)  shows  that  in  the  limit  of  a  vanishingly 
small  crack  length 
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*  =  Jl  -  (a/R)2 


(3.11) 


This  is  a  highly  interesting  result  which  though  precisely  valid  only  for  small  values  of 
"a",  it  nonetheless  satisfies,  for  circular  regions,  the  condition  0  =  0  for  a  =  R. 
Thus,  when  the  crack  has  spread  across  the  entire  circular  region  —  0  as 
expected.  This  suggests  that  in  a  strip  of  finite  width  b  and  height  h 

=  f(l  -  (a/b)2,  h/b)  (3.12) 


where  the  form  of  the  function  f  is  to  be  determined  from  an  exact  elastic  solution  of 
the  problem. 


We  now  return  to  the  Griffith  crack  criterion 
entire  circular  region  as  a  thermodynamic  system,  f 
system  and  since  ^  depends  on  "a",  then  so  is 
generalized  strain  e  where 


i.e.,  Eq.  (3.6).  Regarding  the 
is  an  internal  variable  of  the 
a".  We  proceed  to  define  a 


e  a  eR/¥ 


(3.13) 


and  a  generalized  stress  a  where 

I  a  =  aRlr  (3.14) 

then 

a  =  a ¥/9e  =  Ee^2  (3.15) 

recovering  the  relation 

a  =  Ee^2  (3.16) 

It  also  follows  from  the  usual  thermodynamic  considerations  that  -  Bt/aa  is  the 
internal  force  operating  on  the  fracture  mechanism,  i.e.,  the  "driving  thermodynamic 
force"  that  is  the  cause  of  crack  growth.  In  view  of  Eqs.  (3.10),  (3.11),  (3.13) 
(3.14).  and  (3.15): 


i 
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(3.17) 


a* 


a* 


9a  [l  -  (a/R)2]2  E 


Now  if  we  use  linear  thermodynamic;  of  internal  variables,  the  classical  equation  of 
evolution  for  "a”  is 


a*/9a  +  T}  da/d£  =  0 


(3.18) 


where  rj  is  the  resistance  to  cracking  and  d£  is  an  intrinsic  time  measure  and  rj  is 
the  resistance  coefficient  of  the  internal  fracture  mechanism.  In  our  case  £  is  the 
damage  coordinate.  If,  therefore,  we  set 


d£  =  da 

Eq.  (3.18)  becomes 
-  a$/3a  =  rj 


(3.19) 


(3.20) 


Thus  Eqs  (3.6)  and  (3.18)  are  completely  reconciled  if  rj  =  27.  Therefore,  using  Eq. 
(3.20) 


a 

c 


iff1  ll  -  (a/R)2! 


(3.21) 


which  in  the  limit  of  vanishingly  small  "a"  agrees  with  the  Griffith  crack  criterion,  but 
is  also  qualitatively  true  for  finite  regions  of  width  R  in  the  range  0  <  a  <  1.  as  will 
be  shown. 

We  have  thus  demonstrated  that  the  griffith  fracture  theory  is  totally  consistent 
with  the  damage  theory  proposed  here,  in  the  context  of  irreversible  thermodynamics 
put  forth  above,  for  the  choice  of  damage  co-ordinate  given  in  Eq.  (3.19).  It  also 
provides  a  means  for  determining  the  dependence  of  ^  on  the  crack  extent  in  Mode  I, 
for  crack  lengths  that  are  small  in  relation  to  the  size  of  the  size  of  the  overall 
thermodynamic  system. 
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4.  GENERALIZATION  TO  OTHER  DOMAINS. 

We  wish  to  extend  the  above  ideas  to  other  domains.  To  this  end  we  begin 
with  a  strip  of  finite  width  b  and  height  h.  Our  purpose  here  is  to  use  linear 
fracture  mechanics  as  a  basis  from  which  to  view  the  damage  theory.  Thus  we  are 
looking  for  a  paradigm,  an  example  of  convergence  of  the  two  theories,  when  the 
material  is  an  elastic  solid  that  fractures  in  an  ideally  brittle  manner. 


In  accordance  with  Eq.  (3.9)  the  free  energy  of  the  plate  (width  2b,  height  2h) 
is  given  by  the  expression 


where 


» 


*  =  l 


2 

1 


(4.2) 


being  the  first  principal  component  of  the  integrity  tensor  Again  we  define  a 
generalized  stress  a  and  strain  e  such  that 

a  =  fffbh  ,  e  =  eJbh  (4.2a,b) 

Thus  the  thermodynamic  relation 

a  =  9$/9e  (4.4) 

gives  rise  to  the  expected  constitutive  relation: 

a  =  Ee$  (4.5) 

Equation  (4.5)  provides  the  first  limiting  condition  on  $  in  the  sense  that 

I im  $  =  1  (4.6) 

h*» 

since  in  the  limit  of  h*»,  a  =  Ee. 

From  the  theory  of  fracture  mechanics: 

-  9*/9a  =  Kj/E  (4.7) 
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This,  in  view  of  Eq.  (3.18)  damage  theory  and  fracture  mechanics  are  brought  into 
correspondence  if 


(4.8) 


and  da  =  d£  in  accordance  with  Eq.  (3.19). 

Using  Eqs.  (4.7)  and  (4.1): 

Ee2(a.b2]i  =  K2/E  (4.9) 

a  2  2 

where  $  =  6$/9(l  -  a  /b  )  at  constant  h/b.  At  this  point,  using  Eq.  (4.5)  we  find 
that. 


In  the  notation  of  fracture  mechanics  (11] 
KI  =  rrJTi  F(|  ,  £) 

Thus 


(4.10) 


(4 . 10a) 


(4.10b) 


Since  F  is  known  from  an  elastic  analysis,  we  have,  therefore,  established  an  exact 
connection  between  the  damage  tensor  and  the  crack  length  "a"  under  Mode  I 
conditions  and  the  stipulation  that  the  damage  co-ordinate  £  is  a  linear  function  of 
the  crack  length. 

2  2 

More  specifically  if  we  write  1  -  a  /b  =  x,  h/b  =  y,  Eq.  (4.10b)  becomes 
'  (*]  =  y  p2(x*y)  (4.10c) 

Thus 

*  =  - ^b-  Q—  (4 .  lOd) 

c  -  tJF^dx 
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where  c  is  a  constant  of  integration  to  be  determined  from  the  condition 


=  0 


(4 . lOe) 


Equation  (4.9)  provides  another  limiting  condition  on  $.  Fracture  mechanics  requires 
that  [11],  in  the  limit  of  h+® 


K2  =  a2*aF2(b/a) 


(4.11) 


and  since  under  these  conditions  a  =  Ee.  it  follows  from  Eq.  (3.9)  that 


(4.12) 


where  the  function  0*  is  a  function  of  (a/b)  only. 

Now  if  we  use  limiting  conditions  (4.6)  and  (4.12),  and  in  view  of  Eq.  (4.10), 
we  find  that: 

Kj  =  ofra  IF  (4.13) 


thus  providing  a  relation  between  the  damage  function  and  the  fracture  toughness 
coefficient  in  Mode  I.  for  a  strip  of  finite  width.  Note  if  a^  is  the  effective  area 
representation  of  the  damage  on  a  plane  normal  to  N  then  <p  (a  )  must  satisfy  the 
end  conditions  0(0)  =  1;  ^(1)  =  0. 

A  Specific  Example. 

To  illustrate  the  ideas  involved  in  the  section  we  consider  a  form  of  $  given  in 
Eq.  (4.14),  whish  satisfies  the  limiting  conditions  $  =  1  at  a  =  0  and  $  =  1  at 
h  =  ».  Thus  in  the  limit  of  h*®  we  set: 


0  = 


1 


rb 

(§)? 


(4.14) 


and  explore  the  extent  to  which  K(  derived  from  this  form,  agrees  with  the  one 
derived  from  the  principles  of  fracture  mechanics. 


Use  of  Eq.  (4.10b)  gives  the  result 


Ffa  k)  =  - _ 1 _ 

r  lb  '  hJ  fa]  2  1/2  (1-nrb/h) 

.  "  IbJ 


In  the  specific  case  of  the  infinite  strip  (h  =  *) 


1 

L 

1  - 

ra1 

Lb. 

CM 

1/2 

(4.15) 


(4.16) 


In  Table  1  we  compare  the  exact  solution  for  F  as  given  in  [10]  with  that  of 
Eq.  (4.16).  One  can  see  that  the  agreement  is  reasonably  close  for  a  trial  solution 
such  as  the  one  given  by  Eq.  (4.14) 


Table  1.  Comparison  of  the  Exact  Values  of  F  with  Those  of  Eq.  (4.16) 


a/b 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

l/F(a/b) 

1 

.994 

.976 

.945 

.901 

.843 

.767 

.672 

.551 

.388 

0 

IMa/b)2]1/2 

1 

.995 

.980 

.954 

.916 

.866 

.800 

.714 

.599 

.435 

0 

Equation  (4.14),  however,  does  not  lead  to  a  good  representation  of  F  for  finite 
values  of  h.  Other  simple  forms  of  $  will  be  explored  in  the  future. 
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5.  MORE  ON  THE  RATE  EQUATION  (3.18). 


In  the  previous  sections  we  showed  a  precise  correspondence  between  the 
integrity  tensor  and  the  stress  intensity  factor  in  the  context  of  linear  fracture 
mechanics.  The  key  to  this  correspondence  is  the  evolution  equation  (3.18)  in  the 
specific  instance  when  the  fracture  coordinate  £  is  such  that 

d£  =  da  (5.1) 

One  disadvantage  of  linear  fracture  mechanics  is  that  the  evolution  law  for  the 
crack  length  degenerates  into  an  instability  statement  in  the  sense  that  when  the 
stress  is  below  a  critical  stress  level  the  rate  of  growth  of  zero,  whereas  at  the 
critical  stress  level  the  crack  grows  in  an  unstable-  fashion. 

Previously,  we  proposed  an  evolution  law  of  damage  given  by  Eq.  (2.2).  where 
the  damage  coordinate  f  is  such  that  when  d£  >  0,  then 

d(  =  kde  ,  de  >  0  (5.2) 

It  is  therefore  of  interest  to  explore  the  effect  of  Eq.  (5.2)  on  the  evolution 
Eq.  (3.18).  In  this  specific  case  we  find  that  the  Griffith  crack  theory  is  used,  in  the 
limit  of  R-»«.  but  in  the  context  of  Eq.  (5.2): 

aire2E  =  (»7/k)(da/de)  (5.3) 

This  equation  may  be  integrated  to  give  the  relation 


kE*  3 


(5.4) 


where  aQ  is  the  initial  crack  length.  Note  that  in  this  approach  that  a  crack  will  not 
grow  unless  an  initial  crack  length  is  assumed.  While  this  result  is  physically 
unrealistic,  (in  the  sense  that  Eq.  (5.4)  does  not  account  for  crack  initiation),  it 
agrees  with  linear  fracture  mechanics  precepts  whereby  a  specimen  with  a  zero  crack 
length  cannot  fail.  Note,  however,  the  rapid  growth  of  "a”  with  e3. 


Other  Forms  of  d( 


The  question  of  coupling  between  "a”  and  e,  just  inevitably  arise  in  so  far  as 
the  definition  of  the  damage  coordinate  £  is  concerned.  In  keeping  with  our  previous 
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ideas  on  the  concept  of  intrinsic  time  we  set  d£  equal  to  the  righthand  side  of 
Eq.  (5.5) 


d£  =  (k2de2  ♦  da2)1/2 


(5.5) 


where  k  is  a  material  constant  that  has  the  dimension  of  length.  Thus  in  view  of 
Eq.  (3.17),  (3.18)  and  (5.5)  and  in  the  limit  of  R+06,  one  finds  that 


da  _ kEfae* _ 

d£  =  Ir,2  -  EWe4)1'2 

Evidently,  unstable  crack  growth  occurs  when  da/de  =  »,  i.e.,  where 
2  2 

7  =  Eare  =  arc  /E 


or  at 


cr 


(5.6) 


(5.8) 


thus  recovering  the  Griffith  crack  condition.  However,  Eq.  (5.6)  is  now  a  crack 
growth  equation  which  describes  the  rate  of  crack  growth  prior  to  instability.  In  view 
of  Eq.  (5.6),  for  low  values  of  e 

da/de  ~  kEirae2/’?  (5.9) 

and  thus  Eq.  (5.4)  applies.  The  crack  grows  continuously  until  fracture  takes  place, 
when  the  denominator  of  the  righthand  side  of  Eq.  (5.6)  becomes  zero,  in  accordance 
with  the  Griffith  condition.  Thus  the  evolution  Eq.  (3.18)  in  conjunction  with  the 
definition  of  damage  coordinate  given  by  Eq.  (5.5)  gives  an  ultimate  Griffith  crack 
conditions,  but  also,  a  crack  growth  law  prior  to  instability. 

These  ideas  will  prove  useful  in  the  future  for  developing  a  comprehensive 
damage  theory  for  brittle  or  semi-brittle  solids. 


37 


> 


6.  EFFECT  OF  SIZE  ON  SPECIMEN  STRENGTH. 

In  Eq.  (3.18)  the  driving  force  on  the  fracture  process  is  the  rate  of  change  of 
the  entire  energy  of  the  specimen  with  respect  to  the  crack  length  (actually  crack  area 
of  unit  width).  In  the  case  of  a  strip,  height  h  and  width  b.  the  fracture  stress  is 
given  by  Eq.  (4.10a).  Thus,  given  two  different  geometrically  similar  specimens  so 
that  a/b  and  h/b  are  the  same  in  both  specimens,  but  different  in  size,  the  larger 
specimen  will  have  a  smaller  fracture  stress  a  (K(  being  a  material  property),  since 

ac  ~  1-fa  (5.10) 

all  other  factors  being  equal.  This  rule  is  not  satisfied  precisely  since  the  inverse 
square  law  is  satisfied  only  approximately  and  only  so  for  larger  values  of  a. 

The  question  of  size  on  fracture  strength  is  basically  an  unresolved  problem 
which  could  be  investigated  in  the  future.  The  effect  can  be  accounted  for,  however, 
in  a  formalistic  manner  by  stipulating  that  rj  (the  damage  resistance  coefficient)  is  a 
function  of  the  geometry  in  some  sense  which  is  not  quantized  at  this  stage.  Thus, 
we  denote  the  geometry  by  G  and  we  proceed  to  set 

V  =  17(G)  (5.11) 

This  formalism  allows  one  to  treat  the  rate  equation  in  terms  of  the  average  density 
f  of  the  free  energy  of  the  specimen. 
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A  THEORY  OF  DAMAGE  IN  CONCRETE 


ABSTRACT 


A  constitutive  theory  of  damage  in  plane  concrete  was  developed  leading  to  two 
different  constitutive  formulations  —  Theories  I  and  II  —  depending  on  the  underlying 
specific  physical  and  geometric  hypotheses.  Theory  II  was  used  to  depict  the  behavior 
of  concrete  in  the  course  of  axial  compression.  Calculated  results  in  terms  of  the 
axial  stress-strain  relation  were  compared  with  their  experimental  counterparts.  Close 
agreement  was  demonstrated  between  theory  and  experiment. 

The  study  yielded  insights  into  the  physics  of  degradation  with  strong  evidence 

that: 

(a)  Substantial  damage  does  not  occur  until  the  material  has  reached 
its  fully  plastic  state. 

(b)  The  rate  of  degradation  is  a  strong  increasing  function  of  the  initial 
elastic  modulus.  This  phenomenon  becomes  particularly  pronounced 
in  the  softening  region  of  the  axial  stress-strain  curve. 

Finding  (b)  has  ramifications  in  regard  to  the  design  of  concrete  structures 
intended  to  survive  damage.  The  inference  is  that  high  strength  concrete  will  degrade 
much  faster  in  the  softening  region  following  the  onset  of  maximum  stress. 

A  earlier  damage  evolution  equation  (Valanis,  1987)  proved  very  effective  in 
describing  the  experimental  findings  in  plane  concrete. 


1.  INTRODUCTION. 


The  effects  of  damage  on  the  constitutive  behavior  of  concrete  pose  difficult 
questions  that  cannot  be  treated  in  a  trivial  manner.  In  the  first  place,  the  physics 
of  damage  in  this  material  is  very  much  a  function  of  the  prevailing  hydrostatic 
stress.  At  low  pressures  a  substantial  amount  of  microcracking  occurs  in  a  brittle 
manner  in  preferred  directions  leading  to  induced  material  anisotropy. 

Higher  pressures  inhibit  brittle  microcracking  and  the  macroscopic  deformation  is 
the  result  of  material  "flow"  in  a  ductile  manner.  At  these  high  pressures  isotropic 
damage  is  present  but  experimental  techniques  to  separate  damage  from  purely  plastic 
constitutive  behavior  are  not  always  explicit.  In  fact  purely  constitutive  theories  of 
concrete  either  ignore  the  presence  of  isotropic  damage  at  these  pressures,  or,  by 
implication,  treat  its  effects  in  a  manner  pertaining  to  plastic  deformation. 

The  nature  of  damage  notwithstanding,  the  underlying  constitutive  behavior  of 
concrete  is  very  complex  in  the  sense  that  history  dependence  and  strain  path  effects 
are  significant. 

Such  observations  are  reported  by  Scavuzzo  et  at.  (1983),  on  the  basis  of 
extensive  experiments  on  medium  strength  plain  concrete.  Specifically  Scavuzzo 
recorded  the  strain  response  to  various  complex  non-proportional  stress  paths  that 
brought  into  focus  the  underlying  influences  of  the  stress  history. 

Recently  Valanis  and  Read  (1986)  used  the  framework  of  the  endochronic  theory 
to  develop  an  advanced  non-linear  constitutive  equation  for  concrete,  that  describes 
correctly  the  strain  response  to  the  paths  investigated  by  Scavuzzo,  in  the  absence  of 
damage. 

In  this  work  constitutive  questions  concerning  the  representation  of  damage  and 
the  incorporation  of  its  effects  in  the  constitutive  equation  of  concrete  are  addressed. 
In  this  context,  the  equation  of  evolution  of  damage  is  of  central  importance.  We 
point  out  at  this  juncture  that  a  constitutive  theory  of  concrete  in  the  presence  of 
damage  (which  in  this  work  we  call  Theory  I)  has  been  applied  to  concrete  under 
conditions  of  axial  tension  (Valanis  and  Read,  1989).  This  theory  proved  too  complex 
for  application  to  three-dimensional  problems.  In  this  work,  therefore,  we  have 
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developed  an  alternative  approach  which  is  much  more  physically  transparent  and  more 
easily  applicable  to  three-dimensional  geometries. 


Damage  in  Concrete. 

In  the  theoretical  treatment  that  follows,  we  maintain  the  position  in  our 
previous  work,  Valanis  (1988,  1989a, b)  to  the  effect  that  the  representation  of  damage 
and  its  evolution  are  two  separate  issues.  The  former  is  of  geometric  character  and 
is  independent  of  the  constitutive  properties  of  the  material  at  hand.  The  evolution 
of  damage  on  the  other  hand,  is  very  much  a  function  of  material  physics  and  is, 
therefore,  a  constitutive  problem.  This  question  is  taken  up  in  greater  detail  in 
Section  2. 

A  substantial  part  of  the  damage  in  concrete  at  low  hydrostatic  stress  occurs  in 
a  "brittle"  manner.  The  work  "brittle”  is  not  used  idly  but  is  meant  to  signify 
damage  consisting  of  microcracks  caused  by  increasing  tensile  strain.  Damage 
increments  take  place  in  preferential  directions  dictated  by  the  eigen-directions  of  the 
tensile  principal  components  of  the  strain  increment  tensor.  Thus  the  directions  of 
microcracking  are  dictated  by  the  eigen-vectors  of  the  tensiie  principal  components  of 
the  strain  increment  tensor. 

The  proposition  regarding  directionality  is  consistent  with  the  observations  that: 
(a)  damage  in  tension  is  in  terms  of  cracks  whose  plane  is  normal  to  the  direction  of 

the  axial  stress:  (b)  in  compression  crack  planes  are  tangential  to  the  direction  of 

the  axial  stress:  (c)  oblique  cracks  in  a  tensile  stress  field  extend  in  planes  normal 
to  the  direction  of  the  axial  stress.  This  we  have  found  in  our  work  of  the  extension 

of  oblique  cracks  in  grey  cast  iron.  See  Part  I  of  this  Final  Report. 

The  above  discussion  on  the  physics  of  damage  in  concrete  will  serve  as  the 
basis  for  the  development  of  a  constitutive  theory  of  damage  in  conjunction  with  a 
constitutive  relation  which  accounts  explicitly  for  the  effect  of  damage  on  the 
mechanical  response  of  this  material. 
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2.  A  THEORETICAL  FRAMEWORK  FOR  DIRECTED-CUM-ISOTROPIC 
DAMAGE  IN  CONCRETE 


THEORY  I. 

The  conceptual  foundations  of  the  present  theory  of  damage  were  laid  in  an 
earlier  work  by  Valanis  as  a  separate  contribution  in  the  work  of  Hegemier,  Read, 
Valanis  and  Murakami  (1987)  (see  Appendix  A).  As  a  result  of  that  work  the 
problem  was  resolved  into  two  separate  parts  one  of  damage  representation  and  one 
of  damage  evolution.  The  former  is  a  geometric  problem  while  the  latter  involves  the 
physics  of  damage  and  is  influenced  by  the  material  constitution. 

In  multiphase  media  such  as  concrete  damage  in  the  form  of  directed  cracks  is 
accompanied  by  the  presence  of  small  randomly  oriented  microcracks  which,  in  the 
mean,  constitute  isotropic  damage  (Diamond  and  Bentur,  1984).  In  Reference  1 
directed  damage  was  represented  by  an  integrity  tensor  value  of  field,  $(x.),  symmetric 
and  of  second  order.  When  $  is  null,  at  some  x',  the  material  has  fully  failed  in  the 
neighborhood  of  x.  and  cannot  support  stress  in  any  direction,  while  if  $  is  unity, 
then  the  corresponding  damage  is  zero.  Other  values  of  J  represent  intermediate 
levels  of  damage. 


To  account  for  the  simultaneous  presence  of  isotropic  and  directed  damage  a 
tensor  field  g(x.)  is  introduced,  which  is  a  measure  of  directed  damage,  and  a  scalar 
field  w(x.),  which  represents  isotropic  damage.  Thus, 


i  = 

(2.1) 

subject  to  the  conditions: 

$(M)  =  t 

(2.2a) 

UO.t)  =  2 

(2.2b) 

To  satisfy  conditions  (2.2a,b)  a  representation  of  the  multiplicative  type  is  adopted 
whereby 

4  =  wg  ;  w(o)  =  1  (2.3) 
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The  effect  of  damage  on  the  constitutive  behavior  is  introduced  by  appealing  to  the 
thermodynamics  of  internal  variables.  In  the  manner  adopted  in  Part  I  of  this  work 
on  damage  in  brittle,  elastic  solids,  the  free  energy  density  of  concrete  is  a  function 
of  $.  However,  because  concrete  is  not  elastic  the  problem  has  to  be  treated 
differently. 

It  has  been  shown  recently  (Valanis,  1989c)  that  the  elastoplastic  process  is  the 
result  of  two  series  processes,  plastic  and  elastic.  Each  process  is  accompanied  by  a 
change  in  the  free  energy  density  (relative  to  the  reference  configuration)  —  f  for  the 
plastic  process  and  f  for  the  elastic  counterpart.  The  total  change  inPthe  free 
energy,  f,  is  the  sum  of  the  two  so  that 

*  =  *e  +  *p  (2-4) 


The  essence  of  a  series  process  is  that  the  same  stress  a  is  the  causal  agent  for 
both  processes  while  the  corresponding  strains  are  additive.  Thus  it  follows  from 
thermodynamics  that 
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ace  a£p 


(2.5) 


while 


(2.6) 


When  the  reference  state  is  a  "natural"  state  devoid  of  initial  internal  stresses  then, 
in  small  deformation  theory,  the  free  energy  is  a  homogeneous  quadratic  function  of 
its  arguments.  Therefore, 


4  =  I  A° 


2  i jk£.  i j  k£. 


(2.7) 


(2.8) 
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The  effect  of  damage  in  the  constitutive  response  is  introduced  by  letting  the 
"elastic  moduli"  Ar,  (r  =  0,  ...  n),  be  quadratic  functions  of  $.  In  the  event  that  in 
the  undamaged  state  &r  are  constant  isotropic  tensors  depicting  a  material  which  is 
initially  isotropic  then  it  is  shown  (Valanis,  1987)  that 


Aijk<>  =  Al*ij*ii  *  A2*ikV 


(2.9) 


where  and  Al  are  the  elastic  constants  in  the  undamaged  state.  In  fact,  since  in 
that  state  $  =  g,  then  the  undamaged  form  of  is  given  by  Eq.  (2.9a) 


Ar 


ijk£  =  Al5ij6k£ 


A25ikV 


(2.9a) 


Damage  During  the  Plastic  Process. 

It  follows  from  Eq.  (2.5)  that 


a.  . 

'J 


(2-10) 


Thus,  upon  substitution  of  Eqs.  (2.3)  and  (2.9)  in  Eq.  (2.10)  one  finds  that 


*  Vik^j  l  [£kl  '  qkJ 


(2.11) 


Though  $  is  positive  semidefinite  in  the  sense  that  g  =  0  in  a  fully  damaged 
material,  for  analytical  purposes  we  will  regard  £  as  positive  definite  in  the  sense  that 
we  can  make  it  approach  zero  to  any  required  degree  of  accuracy  without  making  it 
exactly  zero.  Let 


4  =  t 1/2 


A.  .  = 
•J 


a=l 


a  a  ,1/2 
n. n .  a  ' 

I  I  *  ft 


z;  1  J 


(2-12) 

(2.13) 
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> 


where  n?  are  the  eigenvectors  of  &  and  d  its  eigenvalues,  always  non-negative. 
{ 12  ® 

Hence  always  exist,  for  all  a.  Following  now  the  work  by  Valanis  and  Read 

(1989),  we  introduce  the  transformed  quantities  a.  ep  and  qr  such  that 


=  A*ep»A^  ,  a  =  A»qrA 

M  N  M  m  '  p  rw  rv  rw 

(2.14a) 

a  =  AcrA^ 

rw  /v 

(2.14b) 

where 

> 

-H 

II 

> 

(2.14c) 

Note  that  the  matrix  A  used  by  Valanis  and  Read  (1989)  differs  from  the  present  by 
an  orthogonal  transformation.  However  both  give  rise  to  Eq.  (2.16). 

Use  of  Eqs.  (2.14a,b)  in  Eq.(2.11)  gives  the  following  relation  between  stress 
and  strain  in  the  transformed  state: 


(2.15) 


Also  use  of  Eq.  (2.8)  gives  the  following  expression  for  the  free  energy: 


o  n 

f 

-i  2  . 

2' 

l§  *ir[‘-Up  -  *r) 

]  +  2  A2 

P  r 

£  -8 

(2.16) 


Since  Eq.  (16)  can  be  reached  by  setting  £  =  6  in  Eq.  (2.8)  we  conclude  that  the 
transformed  state  is  one  where  only  isotropic  damage  has  taken  place. 

Equation  (2.16)  provides  a  convenient  vehicle  for  determining  the  material  plastic 
response  in  the  presence  of  directed  damage.  The  evolution  Equation  for  cjr  is 
Eq.  (2.17): 


i 


(2.17) 
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The  resistance  tensors  k,r  must  decrease  with  volumetric  damage  and  must  got  to 
zero  as  w-K).  Thus  we  set 


(2.18) 


where  F  and  F_  are  hydrostatic  and  deviatoric  hardening  functions  respectively. 

n  U 

Solving  Eq.  (2.19)  for  £r  in  the  light  of  Eqs.  (2.20),  (2.21)  and  (2.22)  and 
substituting  in  Eq.  (2.5)  the  following  constitutive  relation  is  obtained  for  a\ 
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(2.23) 


) 


fir 


2  = 


r  de^  r 

^Ezd  -  2  )  dT  *  i 


,1 


Kfzu  - 


dekk 

2']-5Td£ 


;o 


where,  as  shown  previously  (Valanis  and  Read,  1986): 


* 

k(2h)  * 


-a  Zr, 

.  r  r  D 
A2  e 


K  e 
r 


-Vh 


£ 

a  =  ,r 
r  b„ 


K  ■  K  ♦  A3/K  *  b3 


K.  =  Aj  .  1  A2 


(2.23a) 

(2.23b) 

(2.23c) 

(2.23d) 

(2.23e) 


However  recent  work  by  Valanis  and  Peters  (1989)  points  strongly  to  the  possibility 
that 


fd  =  fh  =  F 


(2.24) 


We  shall  adopt  this  view  in  this  work  in  which  case  Eq.  (2.23)  becomes 


“2*  = 
u 


p(*  -  z')dT  dz'  +  2 


K(z  - 


kk 

z')d  dz- 


(2.25) 


where 


d2  =  r 


(2.26) 


At  this  point  we  recast  Eq.  (2.25)  in  the  original  configuration  to  obtain  the  following 
constitutive  relation  in  the  presence  of  directed  and  isotropic  damage: 

9 


I 


,z 

^2  ars  =  AriAsj  ^(z  "  z”)dIT  (AikAj£ekJ 
*  ^ 


^rs  K<2  -  2'>dT  take]*2' 


(2.27) 


Special  Forms  of  Eq,  (2.27). 


Of  particular  interest  is  the  case  where  shear  stress  and  strains  are  absent.  In 
this  event  g.  ep  and  g  (and  A)  are  diagonal  tensors  whose  elements  are  denoted  by 
a.,  e?  and  0.  and  A.)  respectively.  In  this  case  Eq.  (2.25)  takes  the  following  form: 


Vi 

u 


,z 

*\  -  2’>df  (*i  £i)d2' 

J  r\ 


♦  K(z 


d 

z  >dF 


dz‘ 

R=1  r  rJ 


( i  not  summed) 


(2.28) 


The  Case  of  Axial  Stress. 

The  experiments  by  Gopalaratnam  and  Shah  (1985)  and  Valanis  and  Khan 
(1989)  were  performed  under  axial  stress  conditions,  i.e., 

ax  t  0.  a2  =  a3  =  0.  et  t  0  .  e2  =  e2  *  0. 

In  this  instance  Eq.  (2.28)  gives  rise  to  the  following  result: 


Vi  *  #1  E<2  -  2')  dF^dz 


(2.29) 


^oei  = 


‘  2'>dT  ta£i)d2' 


(2.30) 
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where  E(z)  and  u(z)  are  related  to  the  functions  p(z)  and  K(z)  in  the  Laplace 
Transform  plane  by  the  following  equations: 


P(P)  = 


Jl£l 


1  ♦  pi/(p) 


K(p)  = 


3(1  -  2pi/(p) 


(2.31a) 


(2.31b) 


where  a  bar  over  a  quantity  denotes  its  Laplace  Transform.  Note  that  scalar  damage 
has  no  effect  on  the  relation  between  an^  ej- 

We  note  that  Eqs.  (2.29)  and  (2.30)  may  also  be  written  in  the  transformed 
configuration  where  they  have  the  following  form: 


ffj  =  u 


t2~  ~ 


E(z  - 


i/(z  - 


K 

Z-)  jJT  dz 


z'>  d - dz 


(2.32) 


(2.33) 


The  Intrinsic  Time  Scale  z. 

Since  the  purely  plastic  behavior  is  formulated  in  the  undeformed  configuration 
where  damage  effects  are  absent,  the  intrinsic  time  z  should  also  be  defined  in  the 
same  configuration.  In  keeping,  therefore,  with  the  concepts  of  endochronic  plasticity. 


d£2  =  de?.de?.  +  x2 
*  ij  ij 


dekk 


(2.34a) 


where  the  Greek  letter  n  is  a  material  parameter.  Use  may  now  be  made  of  Eq. 
(2.14a)  to  express  dz  in  the  present  configuration.  As  a  result. 


d£2  =  l  .  e?.ep 
3  rrr  sj  ij  rs 


+  n 


*k£.ek2. 


(2.34b) 
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THEORY  II. 


Theory  II  is  motivated  by  a  desire  to  provide  an  alternative  physical  premise  to 
the  mechanism  of  the  plastic  process  and  to  arrive  at  an  analytically  simpler  damage- 
sensitive  constitutive  equation  without  sacrificing  the  rigor  of  derivation  on  the 
physical  properties  of  damage  adopted  in  Theory  I. 


As  previously,  the  point  of  departure  is  the  free  energy  f  given  by  Eq.  (2.8) 
with  the  damage-sensitive  stiffness  tensors  Ar  given  by  Eq.  (2.5).  Equations  (2.10) 
and  (2.11)  follow  so  that  Eqs.  (2.8)  -  (2.11)~are  common  to  both  theories.  However, 
the  resistance  tensors  br  are  not  the  same  in  both  theories.  In  Theory  II  br  are  not 
given  by  Eq.  (2.21,b)  as  other,  additional  conditions,  apply. 

The  basic  difference  lies  in  the  equation  of  evolution  of  the  internal  variables  in 
the  damaged  state.  In  Theory  II  this  has  the  form 


(2.35) 


i.e.,  it  is  expressed  in  the  present  state  rather  than  the  transformed  state.  The 
premise  that  the  effect  of  damage  is  geometric  and  does  not  affect  the  physics  of  the 
plastic  process  is  also  maintained  in  Theory  II.  The  constraint  of  this  premise  on  Eq. 
(2.35)  is  discussed  below. 

In  view  of  Eq.  (2.8),  Eq.  (2.35)  becomes 


-  A 


r 


(2.36) 


Since  's  positive  definite  it  has  a  unique  inverse, 
tensor"  Ar  as  in  Eq.  (2.37),  where 


We  form  the  "relaxation 


(2.37) 


The  underlying  physical  justification  for  the  term:  "relaxation  tensor"  will  become  clear 
in  what  follows. 


In  the  undamaged  isotropic  state  Ar  has  the  form 


Ar.,0  =  AT  <5.  .6,  0  ♦  6 .  ,6 . 0 

i  jkt  1  ijk£.  2  ikj£. 

In  fact,  with  reference  to  Eqs.  (2.23c,d) 
*r  =  A2’ 

Pr  =  3Aj  ♦  Aj 


(2.38) 


(2.39a,b) 


Thus  Aj  and  A^  are  the  determining  constants  of  how  fast  the  material  relaxes,  or 
how  its  memory  fades,  with  respect  to  the  intrinsic  time. 

The  premise  that  the  physics  of  the  plastic  process  is  unaffected  by  damage 

leads  to  the  conclusion  that  the  relaxation  tensors  Ar  remain  constant  during  damage. 

Equation  (2.37)  then  provides  the  key  on  how  the  resistance  tensor  br  are 

affected  by  damage,  since  a  rearrangement  of  this  equation  gives  the  result: 

fer  =  6r  •  [tr]  'l  P-4°> 


Analysis. 

Use  of  Eqs.  (2.36)  and  (2.37)  gives  the  following  evolution  equation  for  the 
internal  variables: 


(2.41) 


At  this  point  we  introduce  the  variable  gr  defined  by  Eq.  (2.42): 


“  9 


r 


(2.42) 


13 


Equations  (2.41)  and  (2.42)  combine  to  give  the  following  evolution  equation  for  gr: 


[,6.  .Pr 

1  i  jrk 


kk 


A2p  • 


(2.43) 


This  equation  is  solved  in  the  normal  fashion  by  taking  its  deviatoric  and  hydrostatic 
parts  in  which  event: 


•J 


\  6.  . 
3  ij 


-pr(z-z')  de 


kk 


dz' 


dz ' 


-a  (z-z')  de. • 

,  '  _ LL 

dz ' 


dz 


(2.44) 


where 


dz 


(2.45) 


£?.  -  i  J.6. . 

ij  3  kk  ij 


(2.46) 


and  F  is  the  hardening  function  of  the  plastic  process. 

At  this  point  we  use  Eqs.  (2.9)  and  (2.10)  in  conjunction  with  Eq.  (2.42)  to 
determine  the  stress  in  terms  of  the  variables  gr.  Hence: 


(2.47) 

where  gr  are  given  by  Eq.  (2.22).  This  completes  the  formulation  of  the  constitutive 
equation  in  the  presence  of  damage. 

A  Simplified  Form  of  Eq.  (2.44). 

The  Poisson  ratio  in  undamaged  materials  --  or  more  correctly,  materials  in  their 
undamaged  state  —  remains  substantially  constant  during  an  axial  test.  Furthermore, 
the  Poisson  ratio  of  concrete  is  small  (about  0.10,  or  less)  as  we  have  found  in  the 
course  of  our  own  axial  co  npression  tests  reported  here  in  a  later  section.  Thus 
smali  deviations  in  the  Poisson’s  ratio  are  not  likely  to  be  of  significance  in  the 
prediction  of  the  constitutive  response  of  concrete. 
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•J 


Alhj  *UL  * 


k2\^K 


An  analysis  of  Eq.  (2.44),  the  details  of  which  are  omitted  here,  shows  that  a 
constant  Poisson  ratio,  requires  as  a  necessary  (but  not  sufficient)  condition: 

fit  =  ar  (2.48) 

and  hence,  in  view  of  Eq.  (2.39a.b) 

Aj  =  0  .  (2.49) 

Equation  (2.49)  reduces  Eq.  (2.44)  to  the  much  simpler  form: 


r 

fi 


r2  -«r(i-o 


J0 


(2.50) 


which,  in  conjunction  with  Eq.  (2.48)  gives  rise  to  the  following  plastic  constitutive 
equation  in  the  presence  of  damage: 


'J 


*.  -4>  - 

(J  K.H 


Aj(z 


<■  \  j  / 

2  )  dT~  dz  + 


ik$i£ 


A9(z 


dek^  .  . 

z  >  5T”  dz 


(2.51) 


where 

A1(z)  =  XZAJe 
a2(z)  =  2ZZA2e  r 


(2.52a) 

(2.52b) 


The  sufficiency  condition  for  a  constant  Poisson  ratio  is  that  Aj(z)  also  be 
proportional  to  A2(z). 

At  this  point  we  may  use  Eq.  (2.3)  to  write  Eq.  (2.51)  in  terms  of  the  isotropic 
damage  u  and  the  directed  damage  Thus: 


2  aij 


de 


Aj(z 


-  O 


k£. 


dz ' 


dz ' 


ik'ie. 


de 


A2(z 


z')- 


kl 


Az 


dz'  (2.53) 
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This  equation  has  an  elegance  and  simplicity  which  is  preferable  to  its  counterpart. 
Eq.  (2.27),  of  Theory  I.  It  also  appears  to  be  easier  to  deal  with  numerically.  Also 
as  in  Theory  I  the  intrinsic  time  scale  is  given  by  Eq.  (2.59) 


dz2  = 


ep  l2 

kJCkXJ 


<f>.  d>. 
r  i  r*J 


deP.dep 


•J 


rs 


(2.54) 


Equation  (2.53)  will  >e  used  in  Section  4  in  the  analysis  of  experimental  data  on 
concrete  in  axial  tension  and  axial  compression. 

The  Constitutive  Equation  (2.53)  Under  Axial  Stress  Conditions. 


A  straightforward  analysis  of  Eq.  (2.53)  under  conditions  of  homogeneous  axial 
stress  (0^  #  0,  0s  =  =  0,  t  0,  t  0)  yields  the  following  result: 


E(z  -  z')  dz' 


*2e2  =  *  ^lel 


(2.55) 


(2.56) 


where 


E  = 


A2  [3A1  -  A2] 

2A,  -  A„ 


(2.57) 


v 


(2.58) 


Evidently,  since  A^  is  proportional  to  A^.  it  follows  that  i/(z)  is  indeed  a  constant 
function. 
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3.  THE  DAMAGE  EVOLUTION  EQUATION. 


The  task  of  determining  the  form  of  the  damage  evolution  equation  in  the 
plastic  state  is  difficult.  If  we  follow  the  approach  taken  in  the  case  of  the  elastic¬ 
fracturing  solid  then  the  driving  thermodynamic  force  that  causes  damage  is  in  the 
case  of  directed  damage  and  Qu  in  the  case  of  isotropic  damage,  where 


(f  =  - 


du 


(3.1a,b) 


noting  that  in  the  case  of  plastic  materials 

*  *  K  *  1*p  (3-2) 

Insofar  as  directed  damage  is  concerned  the  evolution  equation  would  be 


(3.3a) 


while  in  the  case  of  isotropic  damage 

dw  =  k  (Qw)ndn 
o 


(3.3b) 


where  dg  is  an  eigenvalue  of  dg,  is  projected  normally  on  the  eigenvector  nr 
and  £r  is  a  damage  coordinate  in  the  direction  of  nr.  Also  D  is  a  damage  coordinate 
for  the  isotropic  damage  process.  The  parameters  k 
constants. 


k  ,  m  and  n  are  material 


This  approach  leads  to  enormously  complex  evolution  equations  involving  both 
the  elastic  and  plastic  parts  of  the  free  energy.  Furthermore  since  g  depends  on 
both  g  and  ui  Eqs.  (3.3a)  and  (3.3b)  are  explicitly  coupled.  Thus  while  the  above 
approach  is  conceptually  appealing,  it  is  analytically  unwieldy. 

A  Useful  Approximation. 

There  are  certain  physically  motivated  approximations  which  lead  to  substantially 
simplified  forms  of  the  damage  evolution  equations  that  can  be  of  practical  value  in 
the  development  of  a  useful  model.  The  first  stems  from  the  fact  that  the  Poisson 


ratio  of  concrete  in  the  undamaged  state  is  small  (less  than  0.1)  and  its  contribution 
to  the  free  energy  may  be  ignorable.  Perusal  of  Eq.  (2.58)  shows  that  for  this  to  be 
so  the  memory  function  A^z)  must  be  zero. 

Also  an  examination  of  the  experimental  data  by  Gopalaratnam  and  Shah  (1985) 
and  Valanis  and  Khan  (1989)  shows  that  significant  damage  does  not  occur  until 
plastic  effects  approach  a  saturated  state,  i.e.,  at  a  value  of  the  plastic  strain  at 
which  the  axial  stress  would  be  close  to  a  saturated  value  had  damage  not  taken 
place. 


Since  during  the  saturated  undamaged  state  the  stress  remains  constant  it 
follows  that  the  elastic  strain  also  remains  constant.  Thus,  under  uniaxial  conditions 
the  physics  is  such  that  damage  —  to  a  very  good  approximation  --  takes  place  at 
constant  elastic  strain.  Furthermore  a  close  analysis  of  the  experimental  data  of 
Gopalaratnam  and  Shah  shows  that  this  in  fact  is  the  case.  The  obvious  immediate 
generalization  of  this  observation  in  three  dimensions  would  be  that  (a)  the  elastic 
free  energy  ^  remains  constant  during  damage  and  thus  /9g  =  0,  0^  /5 u  =  0, 
and  (b)  the  agent  of  damage  is  essentially  the  plastic  strain  tensor  ep. 

Direction  Damage. 

Addressing  the  case  of  directed  damage  first,  it  follows  that 


=  A„p  .  p  .  <t .  . 

2rrir sjri j 


(3.4) 


The  fact  that  significant  damage  begins  to  occur  at  the  onset  of  the  saturated 
plastic  state,  means  that  gr  would  be  near  saturation  where  the  condition  of 
dgr/dz  =  0  would  be  satisfied  closely,  if  not  exactly. 

The  result  of  that  idealization,  in  conjunction  with  Eq.  (2.48)  and  the  fact  that 
Aj  =  0.  is  that  during  damage, 

A2PU  =  dei>j/dz  (34a) 


Insertion  of  this  equation  in  Eq.  (3.1)  gives  the  relation: 


(3.4b) 


.  dep.  dep. 

-  A  — —  SJ  A 
\s  A2  dz  dz  fij 


The  actual  evolution  equation  is  Eq.  (3.3a),  i.e., 


d*r  *  -  kWr>)'V 


(3.5) 


where 


,<o  . 


•J  '  J 


(3.5a) 


If  we  now  make  use  of  Eq.  (3.3)  in  Eq.  (3.5), 
/  dep  de?.  r  \  r  \ 

«  r)  ■  *2  -t  Kt  "[r)"jr) 


At  this  point  we  note  that  n|r*  are  the  eigenvectors  (orthogonal)  of  dep.  This  fact 
reduces  Eq.  (3.6)  to  the  form: 


(3.6) 


,  »  der  der 

,(r)  P  p  r  r 


dz  dz 


(r  not  summed) 


(3.7) 


At  this  point  we  note  that,  whenever  Eq.  (3.5)  applies,  i.e.,  condition  (2.39)  is 
satisfied. 


and 


der  =  d£r 
P  s 


^  ■  u°ir)4r> 


(3.8) 


(3.9) 


i.e.,  ^  is  the  normal  projection  of  g  in  the  direction  n^.  Hence  if  we  let 


k  =  k 
r 


id 

dz  . 


(3.10) 
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then  evolution  Eq.  (3.4)  becomes,  simply. 


v  *  -  krK;r)  Nr 


(r  not  summed) 


Isotropic  Damage. 


In  a  similar  manner,  since  is  constant  during  damage: 


9  u 


(3.11) 


(3-12) 


Terms  in  Eq.  (3.3b)  are  now  defined  except  for  still  the  explicit  coupling  between 
and  Qu  persists  leading  to  complicated  forms  which  are  difficult  to  handle  at  this 
stage  of  theoretical  and  experimental  development.  For  this  reason  we  develop 
uncoupled  equations  from  an  idealization  to  the  effect  that  £  is  oblivious  of  u  and 
vice-versa,  i.e.,  £  =  5.  Hence  and  Qw  are  now  derived  from  the  relations: 


(3.12a) 


(3.12b) 


The  isotropic  damage  coordinate  fl. 

#  Since  to  a  good  approximation  the  agent  of  damage  (both  g  and  u)  is  the 
plastic  strain  tensor  ep,  it  follows  that  dfl  must  be  an  isotropic  function  of  the 
increment  in  plastic  strain  dep.  The  physically  most  appealing  such  function  is  the 
quadratic  function.  But  such  a  function  is  the  intrinsic  line  dz.  Hence  we  set 

•  dfi2  =  defjdefj  ♦  a2[deP;]2  (3.13) 


Thus  dfl  is  a  plastic  path  increment  in  the  plastic  strain  space  of  the  damage 
material. 

This  form  of  the  theory  of  evolution  of  damage  is  now  complete.  We  summarize 
below  the  relevant  equations. 
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SUMMARY  OF  EQUATIONS 


Directed  Damage. 


*2  =  krWTder 


nr  r  rn<b 

Q,  =  n-n.Qr. 

<t>  '  J  'J 


/der;derp>0,  er>0, 
dCr  =  <  P  P 

\0;  otherwise 


rp  r  r  p 

e  y  -  n.n.  e~. 
n  i  j  U 


(r 


d^r  a  eigenvalues  of  djg 
der  =  eigenvalues  of  dcp 

p 

nT  =  eigenvectors  of  dep 


k  =  k 
r  1 


.dz  . 


>  3 1 


W=1 


not  summed) 


<*J>° 


not  summed) 


Isotrooic  Damage. 


Discussion. 


The  foregoing  work  on  damage  evolution  was  developed  quite  late  in  the  life  of 
the  present  project  and  therefore  could  not  play  a  role  in  the  theoretical  treatment  of 
the  experiments.  A  simpler  theory  of  damage  evolution  developed  previously  by 
Valanis  (1987)  was  actually  used.  We  felt,  however,  that,  for  completeness,  the 
above  theory  is  of  sufficient  moment  and  depth  to  merit  inclusion  in  the  report. 
Moreover,  its  inclusion  invites  two  points  of  discussion: 

(i)  It  demonstrated  that  the  equation  of  evolution  of  damage  is  not  a 
trivial  matter  and  a  significant  amount  of  experimentation  is 
necessary  for  its  full  understanding  —  particularly  in  the  presence  of 
complex  strain  paths. 

(ii)  It  points  to  possible  directions  for  future  research  in  the  theory  of 
evolution  of  damage. 

A  Simpler  Theory  of  Damage  Evolution. 

A  simpler  evolution  equation  for  both  directed  and  scalar  (isotropic)  damage  was 
developed  earlier  by  Valanis  (1987)  and  was  used  in  applications  to  elastic  materials. 
Though  the  results  showed  the  right  qualitative  trends  the  lack  of  experimental  data 
prevented  any  precise  conclusions  as  to  its  suitability.  The  experimental  data  reported 
in  Part  I  of  this  report  showed  that  the  proposed  form  was  inclined  to  overestimate 
damage  in  regions  of  low  strain  intensity.  For  this  reason  a  more  suitable  equation 
was  developed.  The  earlier  work  on  the  simpler  theory  is  appended  to  the  report. 
See  Appendix  A. 

The  case  in  plastic  materials,  however,  is  different  in  the  sense  that  damage  is 
more  widely  distributed  around  areas  of  geometric  discontinuities.  Thus  we  expect 
that  the  simpler  theory  may  be  better  suited  for  concrete  as  opposed  to  brittle  solids 
such  as  cast  iron  or  glass.  Here  we  shall  be  dealing  with  the  simpler  theory  in  the 
uncoupled  form. 

The  Case  of  Directed  Damage. 

The  evolution  equation  is  the  following: 
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4.  THE  CONSTITUTIVE  RESPONSE  UNDER  AXIAL  COMPRESSION. 


The  Lateral  Strain  Response. 

The  lateral  expansion  of  plane  concrete  of  low  strength  (~2  x  10  psi)  was 
measured  in  the  course  of  an  axial  compression  test  as  described  in  APPENDIX  B  on 
the  experimental  program.  Low  strength  concrete  brought  the  tests  within  the 
capacity  of  the  testing  machine  and  obviated  the  danger  of  sudden  and  catastrophic 
splintering  of  the  specimens  associated  with  highly  brittle  fracture.  This  was  of 
concern  to  the  personnel  making  the  lateral  strain  measurements. 

We  begin  with  a  descriptive  account  of  the  nature  of  the  lateral  expansion  of 
concrete  under  axial  compression.  During  the  beginning  stages  of  the  test  when  the 
compressive  stress  was  small,  the  observed  lateral  expansion  was  negligible  and  hardly 
measurable  with  a  displacement  gauge  of  sensitivity  of  1"  x  10’4.  During  this  stage 
of  the  test  damage  was  minimal  (non-observable)  i>  *s  of  visible  microcracks  on 
the  specimen  surface.  The  Poisson  ratio  could  not  be  measured  accurately  at  this 
stage  but  appeared  to  be  of  the  order  of  0.05. 

As  the  axial  stress  was  increased  visible  damage  in  terms  of  small  micro-cracks 
began  to  appear  in  the  specimen  surface  while  simultaneously  the  lateral  strain  began 
to  increase  and  eventually  became  of  the  same  order  as  the  axial  strain.  This 
coincided  with  an  attainment  of  maximum  compressive  stress  on  the  axial  stress 
strain  diagram. 

Further  increase  in  the  compressive  strain  beyond  the  point  of  maximum  stress 
began  to  cause  substantial  damage  in  the  form  of  long  vertical  cracks  running  along 
the  length  of  the  specimen.  This  coincided  with  a  substantial  drop  in  the 
compressive  stress.  This  region  of  strain  is  known  in  the  literature  as  the 
"softening”  region.  As  the  axial  strain  was  increased  further  into  the  softening  region 
the  lateral  strain  increased  dramatically  and  became  large,  of  the  order  of  three  to 
four  times  that  of  the  axial  strain. 

Thus  excessive  lateral  expansion  in  the  course  of  axial  compression 
is  directly  related  to  the  occurrence  of  substantial  of  damage. 

The  test  was  terminated  when  damage  was  so  severe  that,  with  further  damage, 
the  specimen  could  no  longer  be  regarded  as  a  continuous  medium. 


To  proceed  with  the  theoretical  prediction  of  the  observed  lateral  strain  we  first 
note  that  the  elastic  response  in  the  presence  of  damage  is  given  by  Eqs.  (4.1)  and 
(4.2)  which  parallel  equations  (2.55)  and  (2.56)  and  can  be  obtained  from  these  by 
setting  E(z)  =  constant.  Thus 


al  = 


.2  c  e 
lEeel 


(4.1) 


»2£2  = 


i/<f> 


e 

1£1 


(4.2) 


where  u  is  a  constant  in  both  the  elastic  and  plastic  response  as  discussed  previously. 
In  view  of  Eqs.  (2.56)  and  (4.2) 

^2e2  =  “  ^lel  (4-3) 


Because  the  axial  strain  e ^  is  compressive  <f>^  =  1.  Thus  the  relation  of  the  lateral 
to  the  axial  strain  is  given  by  Eq.  (4.4): 

*2e2  *  "  v  C1  C4-4) 

For  the  purpose  of  probing  the  micromechanical  implications  of  the  integrity 
tens  or  £  (in  this  case  ^ )  a  series  of  tests  was  carried  out  whereby  specimens  with  a 
square  cross-section  were  tested  in  axial  compression  —  with  different  initial  plane 
"through"  cracks  (or  slots)  of  rectangular  area  each  parallel  to  one  longitudinal  side  of 
the  specimen.  Three  different  cracks  were  used  each  of  length  of  1  in.,  2  in.,  and 
3  in.,  respectively.  See  Appendix  B  for  details. 

The  question  that  we  posed  is  whether  the  presence  of  these  cracks  constitutes 
damage  in  which  case  ^  would  be  different  from  unity  at  the  beginning  of  the  test. 
If  that  were  the  case  then  the  relation  between  and  e 2  would  be  different  in  each 
case  and  would  be  correlated  with  the  initial  length  of  the  corresponding  crack. 

Figure  4.1  shows  that  no  such  correlation  exists.  One  of  two  possibilities 
presents  itself:  either  the  cracks  do  not  constitute  damage  for  the  stress  history  in 
question,  or  the  material  constants  are  such  that  the  effect  of  the  initial  cracks  is 
concealed.  Of  course  both  the  above  are  possible. 
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AXIAL  STRAINS  0000 


Figure  4.1.  Lateral  response  of  concrete  i 


axial  compression. 


On  the  other  hand  there  is  an  excellent  correlation  between  the  initial  elastic 
modulus  of  the  specimens  and  the  lateral  strain  attained  in  the  course  of  the 
experiment,  assuming  that  the  initial  cracks  have  no  effect.  Specifically  "stiffer” 
specimens  expanded  laterally  faster  than  their  less  stiff  counterparts.  Thus  there  is 
an  excellent  correlation  between  the  value  of  the  initial  elastic  modulus  and  the  extent 
of  the  lateral  expansion  in  the  course  of  damage.  This  observation  is  consistent  with 
the  degree  of  damage  observed  —  in  terms  of  the  loss  of  strength  following  the 
attainment  of  maximum  stress  during  compression.  The  rate  of  softening  ---  i.e., 
degradation  --  was  faster  for  the  stiffer  specimens. 

Thus  we  conclude  that  the  damage  rate  constant  k,  is  a 
monotonically  increasing  function  of  the  initial  elastic  modulus.  In 
fact  we  found  the  relation 

.  ko(sinh/Ee)"-2  (4.5) 

where  k  is  the  damage  constant  in  the  damage  evolution  Eq.  (4.6) 
and  k^,  ft  and  m  are  material  parameters. 

The  Damage  Evolution  Equation. 

For  reasons  which  were  discussed  earlier  in  Section  3,  we  chose  the  damage 
evolution  equation  (4.6)  which  was  originally  proposed  by  Valanis  (1986).  Thus 

d*2  =  -  yV6  2  («•<» 


where  k^  and  m  are  material  parameters.  Since  is  a  monotonically  increasing 
function  of  e  ,  the  parameter  m  must  satisfy  the  inequality  m  >  1  as  can  be  seen 
directly  from  tq.  (4.7)  (see  Valanis  (1986),  Appendix  A).  Integration  of  Eq.  (4.6) 
the  light  of  an  initial  value  <p  gives  the  result: 


in 


We  propose  this  equation  as  the  basis  for  the  prediction  of  the  observed  results 
shown  in  Fig.  6.1.  where  axial  and  lateral  strains  are  shown  in  units  of  10  4.  As 
previously  we  shall  use  10  as  the  unit  of  strain.  We  note  that  since  e2  does  not 
correlate  with  tp 2Q  then  k  must  be  a  large  number,  if  ^2Q  is  of  order  of  unity.  This 
turned  out  to  be  the  case.  Equation  (4.7)  then  reduces  to  the  form: 


_  _L_ 

*2  =  {y”-i>i£2i  l’""1 


(4.8) 


which  is  approximate  but  represents  the  real  situation  closely.  Equation  (4.8)  in 
conjunction  with  Eq.  (4.3)  gives  the  following  expression  relating  e2  to 


1 

ry»  - 1)]""2 


m-1  m-1 


m-2 


(4.9) 


or, 


c  e 


m' 
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where 


c 


1  m-1 


m' 


m-1 

m-2 


(4.10) 

(4.10a) 

(4.10b) 


Equation  (4.10)  was  found  to  describe  the  experimental  data  of  Fig.  4.1  in  excellent 
fashion  and  a  relation  between  c  and  Eg  was  thereby  found  and  shown  in  Figure  4.2. 
The  analytical  representation  of  this  relation  is  given  by  Eq.  (4.11): 

c  =  cq  si nh  p  E@  (4.11) 


where  c  =2.2 
o 


x  10‘6  and  p  =  1.310  x  10-6' 
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Thus 


*  kjsinh  p  Ee)m'2 


(4.12) 


where 


k  =  - 7 — 77-  c  m-2 

0  (m  -  0 


m  =  2.36 


Discussion  And  Conclusions. 


(4.13a) 


(4.13b) 


The  experiments  on  plane  concrete  in  conjunction  with  the  theory  of  damage 
developed  herewith  lead  us  to  the  following  conclusions: 

The  damage  induced  by  vertical  cracks  in  concrete  specimens  in  axial  (vertical) 
compression: 

(i)  DOv'  not  have  an  effect  on  the  axial  stress-strain  relation  in 
accordance  with  the  predictions  of  the  theory. 

(ii)  Initial  damage  in  the  form  of  vertical  cracks,  is  either  negligible 
or  the  material  constants  are  such  that  the  lateral  strain 
response  is  unaffected  by  its  presence. 

(iii)  The  theory  predicts  the  lateral  strain  response  due  to  axial 
compression  in  excellent  fashion  and  in  so  doing  establishes  an 
important  relation  between  the  initial  elastic  modulus  and  the 
rate  of  directed  damage.  Specifically  the  higher  the  initial 
modulus  the  higher  the  rate  of  damage.  It  is  significant  to 
note  that  in  elastic-fracturing  solids  failure  is  immediate  with  an 
ideally  infinite  rate  of  damage. 

The  Axial  Stress-Strain  Response. 

Commonly,  the  axial  compression  test  is  treated  as  a  one-dimensional  problem  in 
terms  of  an  axial  stress-strain  relation.  However  while  the  stress  tensor  has  only  one 
non-vanishing  diagonal  component  the  strain  tensor  has  three  such  non-vani:  ’  ing 
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components.  In  fact  the  lateral  strain  response  is  in  itself  a  difficult  constitutive 
problem  when  damage  in  present  —  as  we  have  seen  in  the  first  part  of  this  section. 


To  treat  the  axial  stress-strain  response  we  recall  Eq.  (2.55).  In  the  light  of 
the  fact  that  the  axial  strain  is  compressive,  f  =  1.  this  equation  becomes 


a,  -  u 


E(z  -  z')  dz 


(4.14) 


Since  in  this  test  the  hydrostatic  pressure  is  small  and  to  keep  analytical 
complications  to  a  minimum,  we  set  F  =  1  in  Eq.  (2.45)  so  that 


►  dz2  =  df2  =  «2(^J2^ir^sd£?j<s  <415> 

In  the  theory  of  endochronic  plasticity  the  value  of  k  is  important  in  shear-volumetric 
interactions  relating  to  densification  or  dilatancy.  Since  no  such  effects  are  present 
here  its  effect  is  not  critical  and  its  value  is  therefore  set  equal  to  zero.  Thus  in  the 
case  of  the  axial  test  Eq.  (4.15)  becomes 


2  2 

dz  ^  d^  =  [dej]  +  2  (^2) 

> 

again  because  ^  =  1. 


(4.16) 


It  may  be  seen  from  Eq.  (4.16)  that  to  calculate  the  axial  stress  response  from 
Eq.  (4.14)  one  must  know  <f>2  and  ei)  as  a  function  of  e^.  To  find  this  functional 
relation  we  recall  Eqs.  (2.56)  and  (4.2)  from  which  it  follows  that 


P 


■_2 

.P 


(4-17) 
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Thus : 


(4.18) 


We  are  interested  in  the  stress  response  in  the  damage  region  in  which  e*  is 
substantially  constant  and  equal  to  its  value  at  the  peak  stress.  Thus  can  be 
found  from  the  measured  strains  shown  in  Figure  4.1. 


We  suspect  however  that  the  term  (^de^)2  is  to  be  small  relative  to  (de£)2  or 
^(de^/dfj)2  to  be  small  relative  to  unity  and  numerical  computations  proved  that  to 
be  the  case.  Thus  for  our  purposes,  during  the  axial  test 


dz  =  d£ 


(4.19) 


which  would  have  been  the  result  had  one  regarded  the  axial  test  as  a  one¬ 
dimensional  problem  —  which  it  is  not.  Equation  (4.14)  now  simplifies  and  becomes 


(4.20) 


where 


M 


E(z)dz' 


J0 


(4.21) 


Thus  to  find  the  experimentally  measured  stress-strain  response  in  compression  as 
depicted  in  Fig.  4.2  one  needs  the  memory  function  E(z)  and  the  scalar  damage 
function  u. 

The  Measured  Axial  Stress-Strain  Response. 

The  measured  stress-strain  response  during  axial  compression  is  shown  in  Figure 
4.2.  Damage  in  the  course  of  the  test  was  minimal  prior  to  attainment  of  maximum 
stress  and  became  visible  in  terms  of  small  vertical  cracks  just  prior 
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to  the  appearance  of  the  peak  stress.  Substantial  damage  in  the  form  of  visible 
vertical  cracks  began  to  appear  following  straining  beyond  the  peak  stress  and  into 
the  softening  region. 

An  important  feature  of  the  stress-strain  curves  of  Figure  4.2  is  the  absence  of 
any  correlation  between  the  initial  and  subsequent  directional  damage  in  terms  of  the 
vertical  cracks,  and  the  measured  stress  response.  This  is  of  course  the  prediction  of 
the  theory  in  accordance  with  Eqs.  (4.10)  and  (4.20).  Thus  in  this  sense,  the  theory 
is  upheld  by  experiment. 

Furthermore,  just  as  in  the  case  of  the  lateral  strain  response,  we  see  that  the 
rate  of  damage  —  in  terms  of  stress  diminution  in  the  damage  region  --  is  an 
increasing  function  of  the  initial  elastic  modulus.  Thus  the  observed  scatter  is  only 
apparent  and  can  be  explained  on  the  basis  of  the  dependence  of  the  rate  of  damage 
on  the  initial  elastic  modulus. 


Observed  Response  in  the  Context  of  the  Damage  Theory. 

To  put  the  observed  data  of  Fig.  4.2  in  an  analytical  context  we  took  the 
following  position: 


(i)  The  plastic  response  of  the  undamaged  state  is  the  same  for  all 
specimens. 

(ii)  The  rate  of  degradation  (damage)  is  an  increasing  function  of  the 
initial  elastic  modulus.  This  was  the  case  in  the  lateral  strain 
response. 


We  note  from  Eq.  (4.20)  that  M(e^)  is  a1(eP)lw_i>  i-e.,  the  stress  response  in 
the  absence  of  scalar  damage.  For  convenience,  we  take  this  to  be  the  extrapolation 
of  the  stress-strain  curve  (prior  to  damage)  of  the  specimen  with  the  maximum  initial 
elastic  modulus.  This  is  shown  as  the  dotted  curve  in  Fig.  4.3a.  We  adopt  the 
endochronic  representation  of  this  curve  by  the  functional  relation 


o 


1 


(4.22) 


where  a  is  close  to  unity.  In  this  case  the  presumed  plastic  stress-strain  curve  of 
the  undamaged  state  is  given  by  Ec|.  (4.22)  with  a  =  0.92  and  =  1.47  ksi,  where 
the  unit  of  strain,  as  before,  is  10  .  The  plastic  response  in  the  undamaged  state  is 
then  calculated  for  the  other  specimens  and  is  shown  by  dottcu  lines  in  same  figure. 


33 


100 


200 


300 


Initial  Elastic  Modulus  (ksi) 


i 

Figure  4.3  Effect  of  initial  elastic  modulus  on  the  directional  and  isotropic  damage 
rate  constants. 
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The  Evolution  Equation  of  Isotropic  Damage. 

Just  as  in  the  case  of  directed  damage  where  we  follow  the  theory  developed  by 
Valanis  (1987).  Thus  the  isotropic  damage  evolution  equation  is  given  by  Eq.  (3.22): 

dw  =  -  k w«n  dfl  (4.23) 

where  fl  is  the  isotropic  damage  coordinate.  We  point  out  a  difference  between  the 

present  notation  and  that  used  previously  in  the  reference  cited  above.  There,  the 

isotropic  damage  coordinate  was  denoted  by  £  and  the  symbol  fl  played  a  different 
role.  We  reason  that  since  u  has  no  directionality  it  is  most  likely  to  be  be 
associated  with  the  plastic  state  of  a  material  in  which  case  dfl  must  be  related  to  a 
"plastic”  strain  increment.  After  a  great  deal  of  search  for  the  right  definition  of  the 
"plastic"  strain  as  the  agent  of  u  we  found  that  the  plastic  strain  of  the  undamaged 
state  was  the  most  appropriate  candidate.  The  physics  of  this  conclusion  is  that  the 
damage  due  to  brittle  fracture  does  not  affect  the  evolution  of  isotropic  damage. 
Thus  we  define  ew  as: 

-<*£-*  <*2U  («-24> 

ISi 


where  de  is  the  "plastic”  strain  increment  which  is  the  cause  of  isotropic  damage,  A 
is  the  elastic  compliance  of  the  undamaged  state,  and  gU  is  obtained  from  Eq.  (2.53) 
by  setting  £  =  6  and  u  =  1.  Furthermore  since  we  are  dealing  with  isotropic 
damage,  it  is  natural  to  set  dfl2  as  a  quadratic  isotropic  function  of  de^..  Thus: 


dfl2  =  a2  [dew  -1 
w  l  il  J 


2  ,  w  ,  w 

+  de . .de  .  . 
•J  'J 


(4-25) 


where  Ay  is  a  material  parameter.  In  the  case  of  simple  compression  which  will  be 
considered  in  Section  4,  we  set 


dew  = 

de-f 

(4.26) 

e 

dfl  = 

ldewl 

(4-27) 

within  an  immaterial  constant  insofar  as  Eq.  (4.27)  is  concerned. 
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Comparison  of  Theory  v;ith  Experiment. 

We  found  that  Eq.  (4.23)  with  n  =  1  was  appropriate  for  describing  the 
isotropic  damage  of  plane  concrete  under  compression.  The  measured  response  shown 
in  Fig.  4.2  was  used  to  determine  experimentally  the  relation  between  u  and  Q  using 
the  definition  of  dfl  given  by  Eqs.  (4.26)  and  (4.27).  This  is  shown  in  Fig.  4.4. 
As  pointed  out  above,  the  relations  between  u  and  fl  (for  various  values  of  the  initial 
elastic  modulus)  satisfy  Eq.  (4.23),  with  n  =  1.  Integration  of  Eq.  (4.23)  in  the  light 
of  the  initial  condition  w(0)  =  1.  gives  the  relation 


n 


(4.28) 


It  is  quite  apparent  that  the  damage  rate  constant  is  a  function  of  the  initial 

modulus  E  and  the  appropriate  functional  relationship  is  shown  in  Figure  4.3 
e 

On  the  basis  of  these  findings  the  theoretically  depicted  response  is  shown  in 
Fig.  4.2.  As  may  be  seen  a  good  agreement  between  theory  and  experiment  has  been 
achieved.  While  the  theoretical  analysis  is  not  totally  predictive  it  does  provide  a 
consistent  conceptual  framework  for  the  interpretation  of  the  experimental  findings. 

Discussion  and  Conclusions. 


A  theoretical  development  in  the  form  of  Theory  II  in  Part  II  of  this  report  led 
to  a  constitutive  formulation  for  concrete  in  the  presence  of  damage  —  in  conjunction 
with  two  evolution  equations  each  motivated  by  two  different  physical  considerations: 
one  relating  to  thermodynamics  and  the  other  to  the  physics  of  annihilation  of 
species.  In  the  analysis  the  second  evolution  equation  was  used  because  of  its 
simplicity  and  the  fact  that  it  lends  itself  more  easily  to  numerical  computations. 

The  constitutive  relation  in  question,  i.e.,  Eq.  (2.53)  and  specifically,  its  one¬ 
dimensional  counterparts.  Eqs.  (2.55)  and  (2.56),  were  put  to  the  test  by  comparison 
with  experimental  data  on  concrete  under  axial  compression.  The  results  of  the 
comparison  of  the  theory  with  experiment  were  the  following: 


(i) 


Both  directional  and  isotropic  damage  are  necessary  to  explain 
the  behavior  of  concrete  under  compression. 


37 


This  is  in  contrast  to  the  case  of  axial  tension  where  isotropic  damage  is 
present,  yet  its  inclusion  in  the  constitutive  behavior  is  not  necessary  to  explain 
experimentally  observed  phenomena  which  can  be  modeled  by  directional  damage  alone 
(see  Valanis  and  Read,  1989). 

(ii)  The  theoretical  prediction,  that  the  stress-strain  response  is 
unaffected  by  directional  damage  but  depends  significantly  on 
isotropic  damage,  was  upheld  by  the  experiments. 

(iii)  The  parallel  theoretical  prediction  that  the  lateral  strain  response  -- 
with  excessive  lateral  expansion  developing  progressively  with  damage  — 
does  not  depend  on  the  isotropic  damage  but  is  governed  by  the 
directional  damage  is  also  upheld  by  experiment. 

(iv)  The  "simpler"  form  of  the  damage  evolution  equation,  associated  with 
species  annihilation,  gives  good  agreement  with  experiment  but  this 
result  is  limited  to  observations  associated  with  macroscopically 
homogeneous  deformation  (and  an  axial  loading  history)  and  further 
experiments  associated  with  strain  gradient  fields  are  necessary  to  verify 
its  range  of  applicability. 

(v)  The  micromechanical  experiments  in  terms  of  initial  vertical  cracks 
designed  to  determine  ^  in  terms  of  crack  geometry  (by  measuring  the 
lateral  response  during  axial  compression)  did  not  help  in  this  instance 
to  establish  such  a  relation  because  the  material  damage  rate  constant 
was  such  that  its  magnitude  concealed  the  effect  of  the  initial  cracks  in 
the  lateral  strain  response.  Further  experiments  in  this  regard  are 
indicated. 

(vi)  The  micromechanical  question  regarding  the  geometric  relation  between 
microdamage  and  $  as  well  as  u  is  still  open.  We  believe  that 
establishing  such  a  relation  will  depend  a  great  deal  on  the  correct 
form  of  the  damage  evolution  equation  which  must  be  corroborated  by 
experiments  with  heterogeneous  strain  fields. 
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APPENDIX  A 


NEW  CONTINUOUS  DAMAGE  MODEL  FOR  CONCRETE 
A.S  INTRODUCTION. 

The  purpose  of  this  appendix  is  to  describe,  in  detail,  a  damage-cracking  model 
that  has  been  developed  under  the  present  program.  The  model  was  formulated  by 
Valanis  (1985)  and  is  presently  under  development.  The  discussion  given  below 
largely  follows  that  given  by  Valanis  (1985)  but  does  not  reflect  various  topics  that 
are  presently  under  study,  which  include  (a)  the  addition  of  plasticity,  (b)  size  effects 
and  (c)  strain  rate  dependence. 

A. 2  FORMULATION  OF  MODEL. 

Consider  a  linearly  elastic  material  which  is  isotropic  in  its  virgin  unstrained 
state  and  undergoes  small  deformation  at  isothermal  conditions.  As  the  deformation 
proceeds,  damage  gradually  develops,  which  in  turn  reduces  the  material  integrity.  For 
sufficiently  large  deformation,  the  accumulated  damage  can  lead  to  complete  fracture, 
in  which  the  material  cannot  support  tensile  stress  in  a  direction  normal  to  the 
damage.  We  shall  formulate  this  damage-cracking  process  within  the  context  of 
irreversible  thermodynamics  and.  for  this  purpose,  we  introduce  the  free  energy  *  per 
unit  volume.  Since  fracture  leads  to  a  reduction  in  material  integrity,  we  introduce  an 
integrity  tensor  jg  which  will  be  discussed  at  length  later  and  which  is  symmetrical 
such  that  t  =  §,  (the  unit  tensor),  when  the  material  is  in  its  virgin  (undamaged) 
state,  and  g  =  0  when  the  material  has  fully  failed,  i.e  .  when  the  material  cannot 
support  stress  in  any  direction.  We  thus  set: 


*  =  *(£,*)  (A-l) 

where  e  denotes  the  strain  tensor.  In  a  thermodynamic  sense,  g  now  plays  the  role 
of  an  internal  variable,  in  which  case  we  can  write  the  following  relations: 


2 


(A-2) 


This  nomenclature  is  preferred  to  free  energy  density,  in  view  of  the  concept  of 
"irreducible"  material  volume  which  must  be  introduced  when  dealing  with 
heterogeneous  multi-phase  materials,  such  as  plain  concrete. 


where  £  is  the  stress  tensor  and  g  denotes  the  internal  force,  dual  to  4  which  is 
driving  the  fracturing  process. 

In  view  of  the  definition  of  the  following  relations  must  hold: 


♦  U,0)  =  0  ,  *(0,4)  =  0  ,  (A-3) 

the  first  meaning  that  a  fully  failed  material  cannot  contain  free  energy  and  the 
second  that  an  unstrained  material  must  have  zero  free  energy:  both  conditions  are 
relative  to  the  reference  state.  Furthermore,  the  following  conditions  must  also  hold: 


) 


e(r.S)  *  s  .  e(8.<)  « 2  (*-«) 

9U.2)  =  2  (*-5) 

Equations  (A-4)  require  that  the  stress  vanish  in  the  fully  failed  material  and  that, 
since  the  material  remains  elastic  in  the  damaged  state,  the  stress  will  vanish  at  zero 
strain.  Equation  (3.5)  stipulates  that  the  internal  fracture  causing  force  g  must 
vanish  in  the  fully  failed  material. 

If  we  now  expand  *  in  a  Taylor  series  in  £.  retain  terms  no  higher  than  the 
quadratic,  and  observe  relations  (3.3).  we  obtain  the  following  form  of  ♦: 


*  = 


Ci  jk£e  i  jCk£. 


(A-6) 


where  C  =  C(4).  It  is  further  assumed  that  C  depends  upon  4  in  a  purely  quadratic 
manner.  Since  the  material  is  assumed  to  be  isotropic  in  its  virgin  state.  C  may  be 
represented  in  terms  of  outer  products  of  the  unit  tensor  g  and  4  Furthermore, 
since  C  is  purely  quadratic  in  4 ■  it  has  no  other  representation  than 


C. 


ijk£.  ^  ^  i  j^k£.  *  ^  ^ik^j£.' 


(A-7) 


in  view  of  the  symmetries  imposed  upon  it  by  the  symmetry  of  £  implied  in  Eq.  (A- 
6).  The  constants  X  and  2 ft  in  Eq.  (A-7)  must  indeed  by  the  Lame  constants  of  the 
virgin  material,  since 


t 


*  This  assumption  has  recently  been  verified  by  Valanis  (1987).  using  a  non¬ 
thermodynamic  approach  involving  a  mapping  from  undamaged  to  damaged 
material. 
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C i j k£.  =  X  S\j6k L  *  2fi6\k6i t 
when  g  =  £.  "Thus,  we  can  write 


(A-8) 


•"I 


*  *  ♦ik 


tjl'ij'kl 


so  that,  in  view  of  Eq.  (A-2a).  g  us  given  by  the  expression: 


(A-9) 


ffij  =  X  ^i  j^k£.ek£.  *  2fi  *ik*j£.Ck £.  (A‘10^ 

This  equation  defines  the  constitutive  response  of  an  elastic  fracturing  material,  once 
we  establish  how  the  integrity  tensor  g  is  related  to  the  fracturing  process.  In 
particular,  an  equation  is  needed  which  describes  the  manner  in  which  g  evolves  with 
accumulated  damage,  and  this  is  considered  below. 

A. 3  EVOLUTION  EQUATION  FOR  ^(DIRECTED  DAMAGE)- 


When  the  material  is  not  strain  rate  sensitive,  an  increment  of  tensile  strain  will 
produce  an  increment  of  damage.  Specifically,  given  an  increment  in  strain.  d£.  let 
dea  be  its  eigenvalues  and  n.a  its  eigenvectors.  If  d£  and  are  coaxial,  then  dea 
will  constitute  an  increase  in  tensile  strain  if  de°  >  0  and  e“>^  0.  where  ea  denote 
the  eigenvalues  of  e.  In  the  more  general  case  where  d£  and  e  are  not  co-axial.  dea 
is  said  to  be  an  increase  in  tensile  strain  if  de°  >  0  and  e.n.n.  =  e  n  >  0  (a  not 
summed).  Now  let  dg  be  the  change  in  the  integrity  tensor  g  due  to  the  increment 
d£.  and  let  n^  be  the  eigenvectors  of  dg.  We  now  specify  that  dg  and  de  are  co¬ 
axial.  so  that  r 


(A -11) 


The  physical  significance  of  Eq.  (A- 11)  is  examined  in  the  discussion  following  Eq. 

(A- 14) . 


•  To  completely  characterize  the  evolution  of  g.  the  relation  between  dga  and  dea 
needs  to  be  specified.  In  this  section,  we  consider  a  type  of  damage  that  is 
completely  ’’directed",  that  is.  it  develops  in  a  specific  direction  without  being 
accompanied  by  non-directed  (or  isotropic)  damage.  This  type  of  damage  is  the 
variety  dealt  with  in  linear  elastic  fracture  mechanics  (LEFM)  and  appears  to  typify 

#  that  which  occurs  in  single  phase  brittle  materials.  Later  on  in  Section  A. 7.  we  shall 
consider  the  more  general  case  in  which  the  advancing  damage  consists  of  both 


A- 3 


1 


•  directed  and  non-directed  (isotropic)  components.  Such  damage  is  characteristic  of 

multi-phase  materials,  such  as  plain  concrete,  as  discussed  by  Diamond  and  Bentur 
(1984). 

Restricting  attention  now  to  completely  "directed*  damage,  we  adopt  the 
®  following  expression  relating  dftt  and  de  : 


-  -  ft )  V 


(A-12) 


where 


kdea  ,  for  dcfl  >  0  and  ca  >  0 

n 


l  o 


otherwise 


(A-13) 


Here,  k  and  m  are  positive  constants,  and  =  ^..n?na.  (a  not  summed).  The  form 
of  Eq.  (A-12)  has  its  basis  in  the  physics  of  annihilation.  but  may  be  changed  to 
provide  a  more  realistic  description  of  damage  accumulation  for  a  particular  material. 


In  view  of  Eq.  (A-ll).  we  can  therefore  write 


d^  .  = 
•J 


-  z  ftf*;*, 


(A— 14) 


#  With  this  evolution  equation  for  jg,  the  constitutive  description  of  the  elastic  fracturing 

solid  is  now  complete,  once  the  material  parameters  X.  p,  k  and  m  have  been 
specified. 


Physical  Interpretation  of  Equations 

The  physical  meaning  of  the  above  equations,  and  their  relation  to  the  fracturing 
process,  is  best  understood  by  reference  to  Eq.  (A-10).  which  relates  the  stress  tensor 
to  the  strain  tensor.  Let  I f  be  the  eigenvectors  of  £  and  let  o  and  e  be  referred  to 
a  system  of  coordinates  X  Then,  the  components  of  a  and  t  on  the  X  coordinate 

(X  ^  ^  fl 

system  can  be  written  as: 


rs 


rs 


c..HrHs 
•J  '  J 


c. .NrNs 
•J  '  J 


(A- 15) 
(A-16) 


A-4 


In  this  system  of  coordinates,  g  is  diagonal,  i.e..  =  0  for  r  i  s.  and  Eq.  (A-10) 

takes  the  form: 


a  =  X  i  i  €  *  2u  i  4  e 

r*  ”re"mn  «n  ~  rrm"sn  (l 


Since  £  is  diagonal,  we  can  write: 


(A-17) 


=  X  ^11  ^^rrCrr  +  ^11*11 


11 


*21  .  ;12  =  2/1  *n*22712 


(A— 18) 


^31  *  ^13  =  2/4  ^11^22*13 

Let  us  now  single  out  the  eigenvector  N1  and  examine  what  happens  to  the  stress 
when  the  eigenvalue  ^{—  vanishes.  Note  from  Eqs.  (A-18)  that.  if  =  0.  it 
follows  that  ffn  =  aJ2  =  a  =  0.  that  is.  the  solid  cannot  support  stresses  on  a 
plane  normal  to  N1.  on  which  S  f  =  0.  The  physical  meaning  of  this  is  that 
the  decrease  of  ^  represents  damage  normal  to  the  plane  of  N  and.  as  such,  it  is  a 

measure  of  the  plane  microcracks  that  have  developed  normal  to  N  and/or  the 

increase  in  the  size  of  such  cracks,  so  that  when  —  0.  a  plane  crack  in  the 

accepted  sense  has  formed  across  the  material  element  on  a  plane  with  normal  N  .  so 

that  the  element  cannot  support  tensile  or  shear  stresses  on  that  plane. 

On  the  basis  of  the  above  observations,  the  physical  meaning  of  Eqs.  (A-ll). 
(A-12)  and  (A-13)  now  becomes  dear.  An  increase  in  tensile  strain  defl  on  a  plane 
with  normal  na  causes  planar  damage  on  that  plane,  where  d^fl  is  given  by  Eq. 
(A-12).  in  accordance  with  the  observation  th*‘  "^  anar  microcracks  form  perpendicular 
to  the  direction  of  the  principal  tensile  strain. 

The  constant  k  in  Eq.  (A-13)  is  a  material  parameter  that  reflects  the  fracture 
resistance  or  fracture  toughness  of  the  material.  Eq.  (A-12)  is  representative  of 
processes  involving  annihilation  of  populations  and  catastrophic  processes  in  systems 
where  the  increment  of  annihilation  is  proportional  to  the  state  of  integrity,  which  is 
given  here  by  the  factor  (^a)m.  This  form  is  representative  of  such  systems,  even 
though  any  monotonically  decreasing  function  f(^“).  such  that  f(0)  =  0.  will  likely  do 
as  well.  Finally,  we  note  that  the  damage  tensor  Q  can  be  expressed  in  terms  of  g 
as  follows: 


A- 5 


D. .  = 
•J 


(A-19) 


Y.  fi  -  (^)2]n“n“ 

a  1  1,1 

The  physical  interpretation  of  the  governing  equations  is  now  complete. 

A. 4  APPLICATION  TO  SOME  SIMPLE  CASES  OF  HOMOGENEOUS 
DEFORMATION. 

In  this  section,  analytic  solutions,  based  upon  the  constitutive  model  described 
above,  are  presented  for  several  simple  cases  of  homogeneous  deformation  including 
simple  tension,  simple  compression  and  simple  shear.  These  solutions  provide  valuable 
insight  into  the  characteristics  of  the  model  and  reveal  its  remarkable  predictive 
capabilities,  despite  the  fact  that  it  involves  only  four  material  parameters. 

A. 4.1  Homogeneous  Triaxial  Strain  Fields. 

As  a  preliminary  development,  consider  the  case  in  which  the  strain  field  is 
homogeneous  and  consists  only  of  principal  strains,  i.e..  =  0  for  i  t  j.  Thus,  de 

and  e  are  always  coaxial  and  g  is  coaxial  with  £  by  virtue  of  Eq  (A- 11).  As  a 
result.  I  is  diagonal.  In  view  of  Eq.  (A-10),  £  is  also  diagonal.  The  principal 
stresses  are  given  below  in  terms  of  the  principal  strains  and  the  principal  values  of 

t 

*  MjE  *rlr  *  2«£l 
r 

o2  =  X  +2  E  *rer  ♦  2 ii  f*e2  (A-20) 

o3  =  X  *r€r  -  2 n  *363 

r 

The  appropriate  evolution  equation  for  £  is  obtained  from  Eq.  (A-14)  and  is  given  in 
terms  of  the  principal  components  of  ^  as  follows: 
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d*i  =  *  (*i)m  d*. 


(A— 21) 


where 


* 


fk  d£j  ,  if  d£j  >0  ,  6.  20 
0  ,  otherwise 


We  now  consider  two  special  cases  of  the  above  equations,  namely,  simple  tension 
and  simple  compression. 

A. 4. 2  Simple  Tension. 


Consider  the  case  of  monotonically  increasing  simple  tension.  In  this  case,  a  = 
a  =  0.  and  integration  of  Eqs.  (A-21)  with  e 2  =  c3  and  the  initial  conditions  f>2  = 
f>3  =  1.  shows  that  ^  always.  Hence,  we  find  that 


*2*2  = 


2(X*/») 


ejh*  '  **iei 


(*-22) 


and 


a 


1 


where 


E  _  C2if3X^ 


(A-23) 


(A-24) 


From  Eq.  (A-21a).  the  evolution  equation  for  f  is 
=  0  , 

subject  to  the  initial  condition  ^  (0)  =  1.  Thus,  the  solution  is 
-k£ 

e  ,  if  m  =  1 

=  JL 

[{l  -  (l-m)ke1}1'm  ,  if  (I)  t  1 
Combining  Eqs.  (A-23)  and  (A-26).  we  obtain  the  expression: 


(A-25) 


(A-26) 


; 


(A-27) 


Inspection  of  Eq.  (A-27)  reveals  that,  for  0  <  m  <  1.  ff  is  not  a  monotonic  function 
of  e  but  reaches  a  maximum  and  then  goes  to  zero  for  finite  For  1  <  m  <  3. 

ffj  again  reaches  a  maximum  and  then  decays  monotonically  to  zero  as  *  •.  For 

m  s  3.  monotonically  increases  to  an  asymptotic  constant  value  as  ■»  •.  while 

for  3  <  m  <  •.  <r  increases  monotonically  with  e^.  The  influence  of  m  on  the 

constitutive  response  is  shown  in  Figure  A. la.  Thus,  in  the  range  0  <  m  <  3.  the 
model  exhibits  softening  in  tension  in  agreement  with  experimental  data. 

The  model  also  unloads  elasticity  with  an  effective  modulus  of  E^.  which  is 
diminished  by  the  accumulated  damage.  Thus,  the  model  has  all  of  the 
characteristics  of  a  fracturing  elastic  solid.  It  follows  from  Eq.  (A-22)  that  <  0 
and  hence  f  =  1.  i.e..  no  damage  is  produced  in  the  transverse  direction  due  to  the 
uniaxial  tension  loading. 


A. 4.3  Simple  Compression. 


Consider  now  the  case  of  monotonically  increasing  uniaxial  compression.  Here. 
=  0.  while  Oj  and  are  compressive.  Again.  Eqs.  ( A- 21)  together  with 
=  e3>  and  the  initial  conditions  f  ^  =  1.  lead  to  the  result  that 
always.  Similar  to  the  case  of  simple  tension,  we  find  that 


*2e2  =  '  uhel  (A'28) 

and 

°x  =  EfJcj  (A-29) 

The  damage  evolution  in  this  case  is.  however,  entirely  different  from  that  for  simple 
tension.  Since  both  dCj  and  are  compressive  at  all  times,  we  find  from  Eq.  (A- 

21)  that  d^  =  0  and  thus  ^  =  1  always.  However,  in  view  of  Eq.  (A- 28 ) .  e2  is 
now  tensile  and  hence  (and  ^3)  will  increase  with  the  consequence  that  damage 
will  develop  on  planes  wnich  are  parallel  to  the  axis  of  compression,  in  accordance 
with  experimental  observation;  this  is  the  so-called  axial  splitting  mode  (Horii  and 
Nemat-Nasser.  1985). 
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*1 

(a)  Stress-strain  relation  in  simple  tension 


(b)  Axial  strain  vs.  lateral  strain  in  simple 
compression 


Predicted  responses  for  simple  tension  and  simple  compression,  showing  the 
effect  of  the  parameters  m  on  the  resulting  behavior. 


S(  A  9'l? 


From  Eq.  (A-21).  we  can  write 


d*2  ♦  d€2  =  0 

the  solution  of  which  is 


(A-30) 


{l  -  (l-m)ke2}1_m  ,  m  *  1 


(A-31) 


Thus,  in  view  of  Eq.  (A-29)  and  the  condition  ^  =  1.  the  axial  stress  strain  relation 

is 


ox  *  E  ,  (A-32) 

so  that  the  response  in  the  axial  direction  is  purely  elastic. 

To  find  the  relation  between  the  strains,  we  use  Eqs.  (A-28)  and  (A-31) 
together  with  the  fact  that  ^  =  1.  It  then  follows  that 


-  ^  = 


-ke2 

/  ,  m  =  1 

1 

e2[i  -  (l-m)ke2]1'm  ,  m  r  1 


e2e 


(A- 33) 


Here,  we  may  distinguish  three  different  cases: 
(a)  0  <  m  <  1 


The  strain  is  not  a  monotonic  function  of  It  reaches  a 

absolute  value  and  then  goes  to  zero  for  a  finite  value  of  e^.  Specifically. 


maximum 


lelLax 


=  -  M-l 

k v  In+lJ 


n*l 


(A-34) 


where  n  =  l/(l-m) . 


Again,  the  strain  is  not  a  monotonic  function  of  e ^  It  reaches  a  maximum 
absolute  value  again  given  by  Eq.  (A- 34 )  for  1  <  m  and  then  goes  to  zero  for 
infinite  e^.  For  m  =  1 


^C1  ^max  ”  i/ek  ' 

where  e  denotes  the  base  of  the  natural  logarithm. 


(A- 35) 


(c)  2  <  m  <  • 

In  this  case.  Cj  is  a  monotonic  function  of  and  approaches  the  following 
limiting  condition  for  increasing  lateral  strain: 


I  im  I  I 

|£l|max 


m  =  2 
m  >  2 


The  above  cases  are  depicted  graphically  in  Figure  A.  1(b). 


(A-36) 


A. 4. 4  Collapse  of  a  Block  Under  Axial  Compression. 

It  follows  from  the  above  discussion  that  a  block  consisting  of  elastic  fracturing 
material  will  collapse  under  homogeneous  axial  compression  if  0  <  m  <  2.  since  e 
has  an  upper  bound  given  by  the  expression: 


'  1 1  ma  x  ~  uk 


-m J 


(A-37) 


The  limiting  cases  are 


.tab  •  *  * 1 


1 1  max 


U  ■  ■  *’ 

The  collapse  stress  may  be  calculated  directly  from  Eq.  (A-32).  with  the  result 


(A-38) 


2-m 

hi*.,  =  b  (jb]1”  (*-*>) 

Thus,  the  theory  predicts  the  collapse  of  a  block  under  axial  compression  due  to 
damage  on  planes  parallel  to  the  axis  of  compression.  Again,  this  is  the  axial 
splitting  mode. 


A-ll 


A. 4. 5  Ratio  of  C^'apse  Stresses  in  Tension  and  Compression. 


At  this  point,  it  is  of  interest  to  explore  the  rationality  of  the  above  results  by 
comparing  the  ratio  ac /ol.  where  oc  and  a1  are  collapse  stresses  in  compression  and 
tension,  respectively,  obtained  analytically,  with  that  observed  experimentally  for 
concrete,  even  though  the  end-grip  conditions  in  the  tests  may  not  be  ideal,  as 
assumed  in  solutions  obtained  here.  Restricting  attention  to  the  case  m  =  1.  it 
follows  from  Eqs  (A- 37)  and  (A-39)  that 


(A-40) 


Since  v  **  0.2  for  concrete,  we  have 


\  :  10  ,  (A-41) 

which  is  close  to  experimental  observation  (Raphael.  1984). 

A. 4. 6  Simple  Shear. 

Consider  now  the  case  of  monotonically  increasing  shear  strain  under  conditions 
of  simple  shear  and  small  strain.  The  configuration  of  interest  is  depicted  in  Figure 
A. 2  Inasmuch  as  the  only  non-zero  strain  component  is  e  we  can  write 


=  2  *12C12 


(A-42) 


*.kW=  (*il*j2  *  ^  i  2^j  lJ  C12 


(A— 43) 


Therefore,  in  this  case,  the  general  constitutive  relation 


"ij  *  k  'i/kH'ki  *  ^  'ik'ji'ki 


(A-44) 


reduces  to  the  form 


2 


Therefore,  we  can  write 


°\2  ~  2[X  *12  4  ^(*11*22  *  *12) ) e 12 

(A-46) 

°ll  =  2(X  *  2^)*11*12612 

(A-47) 

°22  ~  2(X  *  2^]*22*12C12 

(A-48) 

For  a  small  shear  strain,  the  maximum  tensile  strain  will  lie  in  a  direction  that 
makes  a  45  degree  angle  with  each  of  the  coordinate  axes.  Consequently,  cracks  will 
develop  in  the  material  as  shown  in  Figure  A. 2.  The  components  of  the  unit  normal 
rj  to  the  crack  are  n^  =  n^  =  l/f2. 

Since  there  is  only  one  cracking  direction,  we  set  a  =  1  in  the  evolutionary 
equation  for  d^..  and  write 


d*ij  *  -  kM"dd  Vj 


(A-49) 


where 


i  =  i. .nn. 
rn  rij  i  j 


de  =  de 


(A-50) 


12 


Since  n^  =  nj  =  1/J"2,  it  follows  from  Eq.  (A-49)  that 
d*ll  =  d*22  =  d*l2 

Hence 


(A— 51) 
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f*12 


dip 


' 1 


11 


d*12  * 


(A— 52) 


since  p.  T.j  initially.  From  Eq.  (A-53)  we  find 

*11  ~  '  *12  ' 

and  Eq.  (A- 50a)  gives 

*n  =  *11  *  *12  ' 


(A-53) 


(A-54) 
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since  =  ^2  *ccording  to  Eq.  (A-51)  . 


Substitution  of  Eq.  (A-53)  into  Eq.  (A-54)  yields 

K  * 1  • 2  *u 

which,  when  combined  with  Eq.  (A-49) ,  leads  to  the  result 

d*ij  *  -  kP  *  2^12]’  df12ninj 

Therefore,  since  n^  =  n^  =  1//2,  we  can  write 
d^12  ~  ~  2  [l  *  2  ^12]  dC12 


(A-55) 


(A-56) 


(A-57) 


which  can  be  integrated  to  give,  for  m  =  1 : 


1  -kf12 

^12  =  2  l*  l. 


(A-58) 


and  for  m  t  1 : 


a  I  ,  1  m-1 

fl2  "  2  k(m-l)e  ♦!  "  1 


(A-59) 


Recalling  Eq.  (A-53)  and  the  fact  that  ^  ^  .  Eq  (A-46)  can  be  placed  in  the 

following  form: 


al2  =  2[(X  4  W12  *  2^12  *  ^]£1 


(A-50) 


Therefore.  Eqs.  (A-58)  to  (A-60)  describe  the  relationship  between  a ^  and  e ^ 
From  Eqs.  (A-58)  and  (A-59)  note  that 


I  i">  ff12  =  ~  2  (X  4  2^)e 


C12*" 


(A— 61 ) 


which  in  view  of  Eq.  (A-60)  leads  to  the  result: 


°12  =  2  *  2^el2 


(A-62) 
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Hence,  when  X  <  2 p.  the  ultimate  slope  given  by  Eq.  (A-62)  will  be  less  than  the 
initial  slope.  2/*.  as  shown  in  Figure  A. 3.  The  case  X  >  2p  on  the  other  hand,  leads 
to  the  ultimate  slope  being  greater  than  the  initial  slope,  i.e..  a  hardening,  which  is 
difficult  to  visualize  from  a  physical  standpoint.  These  limiting  cases,  however,  fall 
outside  of  the  assumption  of  small  strains  made  here  and  therefore  may  not  be 
physically  meaningful. 

A.5  APPLICATION  TO  NON-HOMOGENEOUS  STRAIN  FIELDS. 

In  an  effort  to  explore  the  capability  of  the  continuous  damage  model  to  predict 
the  proper  trends  for  fracture  problems  involving  non-homogeneous  strain  fields,  we 
developed  a  computer  program  for  the  model  and  incorporated  it  into  a  finite  element 
code.  We  then  conducted  a  limited  numerical  study  of  a  fracture  problem,  having  a 
non-homogeneous  strain  field,  for  which  a  theoretical  solution  is  available  from  linear 
elastic  fracture  mechanics  (LEFM).  The  problem,  depicted  in  Figure  A. 4(a).  consists 
of  a  long  flat  plate  of  width  w.  which  has  a  crack  of  length  2a  and  is  acted  upon  by 
a  uniform  tensile  stress  a.  The  problem  is  to  find  the  magnitude  of  a  which  causes 
the  crack  to  grow  catastrophically  and  thereby  cause  fracture  of  the  plate. 

The  finite  element  grid  used  in  the  numerical  study  consisted  of  uniform,  initially 
square  zones  of  unit  length.  Because  of  symmetry,  only  the  right  half  of  the  problem 
was  considered.  Vertically,  there  were  21  zones  while  16  zones  were  used  in  the 
horizontal  direction.  Therefore,  the  width  w  had  the  value  32.  The  effect  of  an 
initial  crack  in  the  x-direction  was  introduced  into  the  grid  by  simple  setting  l  =0 
in  those  zones  through  which  the  crack  passed.  In  the  study.,  three  differentiates 
of  initial  crack  length  were  considered,  namely,  a  =  1.  3  and  6. 

The  results  from  the  numerical  study  are  shown  in  Figure  A. 4(b).  where  the 
stress  af  to  fracture  the  plate  is  plotted  against  the  dimensionless  ratio  a/w.  Also 
shown  on  this  figure  is  the  corresponding  theoretical  prediction  from  LEFM  (see 
Knott.  1973).  which  is  given  by  the  expression: 


where  K)(..  the  mode  I  fracture  toughness,  was  set  to  the  value  7.07  so  that  the 
theoretical  curve  would  pass  through  the  central  numerical  point  shown  in  the  figure. 
As  an  inspection  of  the  figure  reveals,  the  continuous  damage  model  predicts  a  decay 
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(a)  Geometric  configuration  of  flat  plate. 
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Response  of  a  thin  elastic  plate  with  a  crack  to  an  applied 
uniform  tensile  stress. 
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of  of  with  increasing  a/w.  which  is  qualitatively  in  agreement  with  the  theoretical 
result  from  LEFM. 


Exact  agreement  between  the  model  prediction  and  LEFM  is  not  to  be  expected, 
however,  for  several  reasons.  First,  the  model  produces  some  damage  in  every  zone 
of  the  grid  from  the  onset  of  the  tensile  loading,  while  in  LEFM  theory  it  is  assumed 
that  the  damage  remains  confined  to  the  original  crack  until  there  is  catastrophic 
crack  extension.  Secondly,  for  the  material  region  just  ahead  of  the  crack,  the  two 
approaches  lead  to  significantly  different  stress  fields.  In  LEFM.  the  material  ahead  of 
the  crack  remains  linearly  elastic  and  undamaged  until  overall  failure  occurs,  while,  in 
the  numerical  simulation  with  the  model,  the  zones  ahead  of  the  crack  experience 
increasingly  greater  damage  as  a  grows.  Thus,  the  stress  fields  in  the  material  just 
ahead  of  the  crack  are  very  much  different  in  the  two  cases.  In  view  of  these 
considerations,  the  differences  between  the  two  predictions  shown  in  Figure  A. 4  are 
not  surprising.  The  general  trends,  however,  are  similar  and  indicate  that  the  present 
model,  when  used  in  conjunction  with  a  finite  element  code,  leads  to  realistic  results. 
Closer  agreement  between  the  two  predictions  could  probably  be  obtained  by  revising 
the  evolution  equation  for  £. 

A.6  POST-FAILURE  BEHAVIOR  AND  CRACK  CLOSURE. 


A 

It  was  shown  in  Section  A. 3  that  if  an  eigenvalue  f  of  £  vanishes  in  a 
material  element  then  the  element  cannot  support  stress  on  a  plane  normal  to  NA. 
the  latter  being  the  eigenvector  of  g  associated  with  ^A.  In  a  structural  sense  the 
element  has  developed  damage  equivalent  to  a  crack  on  the  plane  normal  to  NA.  It 


follows  that  the  element  has  developed  kinematic  freedom  in  directions  normal  and 
tangential  fo  NA  in  the  following  sense.  Let  Ae be  the  components  of  the  strain 
increments  A£  defined  relative  to  a  system  of  coordinates  xa  coincident  with  the 
eigenvectors  Na  of  Thus 


iTll 

Ae12 

Al13 

Ac22 

Ae23 

4'31 

Ac32 

Ae33. 

(A-64) 


Let  A  *  1.  Perusal  of  Eqs.  (A-18a.b.c)  show  that  if  =  0.  then  not  only 

are  =  c ^  as  a ^  =  ?2J  =  0  as  discussed  previously  but  the  strains  e 

7  =  e  c. .  =  e  become  indeterminate  from  the  constitutive  properties  of  tne 

1*  *1  41  14 

element. 

A.7  DIRECTED  AND  ISOTROPIC  DAMAGE. 

It  has  recently  been  observed  (Diamond  and  Bentur.  1984)  that,  in  multi-phase 
materials,  such  as  plain  concrete,  the  formation  of  directional  cracks  is  usually 
accompanied  by  the  presence  of  small  cracks  randomly  oriented,  which  from  a 
statistical  viewpoint  amount  to  isotropic  damage,  denoted  here  by  Q.  In  this  sense  0 
reduces  the  stress  carrying  capacity  of  a  material  by  the  same  degree  in  all  directions. 
Thus  the  integrity  tensor,  now  denoted  by  must  be  reduced  by  the  presence  of  Q 
and  hence  it  must  be  a  state  function  of  Q.  Thus,  in  symbolic  form: 


In  *  <A-65> 

subject  to  the  following  conditions  which  must  always  apply 

In (o.o)  *  S 

(A-66) 

•  s 


where  £  is  the  unit  matrix. 

To  satisfy  conditions  (A-66)  we  adopt  a  relation  of  the  multiplicative  type,  i.e.,  we 
set 

lfl  •  «(0)|((J  (A-67) 

such  that  «(0)  =  1.  «(1)  =  0.  Furthermore  u  must  be  a  monotonically  decreasing 
function  of  Q.  A  suitable  function  that  satisfies  the  above  constraints  is 

v(D)  =  (1  -  Q)a  (A— 68) 


where  a  i  0.  Thus  the  integrity  tensor  £.  considered  in  the  previous  sections,  is  now 
modified  by  the  presence  of  isotropic  damage  as  specified  by  Eqs.  (A-67)  and  (A-68). 


The  constitutive  equation  thus  becomes 


*  *  2  *2{X[^ijeij)2*  * 

and 

«ij  =  *  2/>  (*-70) 

A.7.1  The  Rate  Equation  for  Q. 

Physical  considerations  suggest  that  isotropic  damage  in  cementitious  materials  is 
caused  mainly  by  shearing  which  gives  rise  to  local  decohesion  (between  the  mortar 
and  the  aggregate,  and/or  the  paste  and  the  sand  particles).  Thus  following  similar 
reasoning  that  gave  rise  to  Eq.  (A-12)  and  noting  that  the  ultimate  value  of  Q  is 
unity,  we  set  dQ  to  be  proportional  to  a  positive  monotonic  function  of  the  isotropic 
integrity  (1-Q).  i.e.. 

dQ  *  f(l  -  Q)df  ,  (A— 71) 

where  dc  is  defined  in  terms  of  the  increment  in  the  deviatoric  strain  tensor  de  .  as 

*  'J 

d£2  =  c2(deijd«ij)  .  •  (A-72) 

and  c  >  0  is  a  material  parameter  A  suitable  form  of  the  function  f  is  as  follows: 

f(l  -  0)  •  (1  -  Q)b  (A-73) 

where  b  >  0.  Solution  of  Eq.  ( A- 71)  in  the  light  of  Eq.  (A-73)  gives  rise  to  the 

expression 


1 

1  -  Q  *  (1  -  (1  -  b)0l"b  , 

which  is  of  the  same  form  as  Eq.  (A-26).  Specifically,  if  b  =  1  then 


(A-74) 


(1  -  Q)  =  e'b?  (A-75) 

We  note  that  c  in  Eq.  (A-72)  plays  the  role  of  an  isotropic  fracture  toughness 
parameter  which  influences  the  rate  of  isotropic  degradation. 


A- 20 


In  the  presence  of  isotropic  damage  the  counterparts  of  Eqs.  (A-20)  are: 


*  **  j 

(A— 76) 

X*2^k‘k  '  *2C2 

(A- 77) 

xr32Irk'k  *  2/i  r3e3 

(A-78) 

In  that  which  follows,  we  illustrate  the  application  of  the  above  equations  to 
two  special  cases  of  interest,  namely,  simple  tension  and  simple  compression. 

(A)  SIMPLE  TENSION 

Here  —  0  and  the  analysis  follows  very  closely  that  of  Section  A. 4  2. 

Specifically,  in  view  of  Eqs.  (A-78): 


^262  ~  ~  v 


(A-79) 


o 


1 


(A-80) 


However,  as  opposed  to  the  analysis  in  Section  A. 4. 2  two  functions  ^  and  u  are  now 
necessary  to  determine  the  axial  stress  a  .  In  the  absence  of  explicit  coupling 
between  directional  and  isotropic  damage,  ^  is  still  given  by  Eq.  (A-26)  in  which 
case 


1 

=  {l  *  (1  -  *")kc1}1"m  (A-81) 

However,  we  expect  that  k  in  Eq.  (A-81)  will  depend  on  fl  since  clearly  the  material 
resistance  to  cracking  will  be  greatly  influenced  by  the  presence  of  the  isotropic 
damage.  Specifically  we  expect  that 

A- 21 


k  =  k(Q) 


(A-82) 


where  k  is  an  increasing  function  of  Q.  In  this  event  Eq  (A-25)  can  no  longer  be 
integrated  directly  and  Eq.  (A-81 )  will  not  hold. 

Evolution  of  Isotropic  Damage. 


In  the  presence  of  axial  loading  (o  =  a.  =  0) 


dS  =  «Jf  d«j  "  d*: 


(A-83) 


which  in  the  light  of  Eq.  (A-79)  becomes 


dS  =  cJI  IdEj  .  *  d[£ 


(*-84) 


where  use  was  made  of  Eq.  (A-12).  Thus,  in  general,  we  have  the  following  coupled 
equations  that  govern  simultaneously  the  evolution  of  directional  and  isotropic  damage 
in  simple  tension: 


d*j  =  -  k(Q)de: 


(A-85) 


dO  =  (1  -  Q)b  cJI  (l  ♦  -  J/k^e1]de1 

Approximate  Solution  of  Eqs.  f A-85)  and  (A-86). 


(A-86) 


To  obtain  a  qualitative  understanding  of  the  simultaneous  effects  of  isotropic  and 
directional  damage  we  solve  Eqs.  (A-85)  and  (A-86)  approximately.  Specifically  we 
uncouple  the  equations  by  setting  k(fl)  equal  to  a  constant  k  and  setting  equal  to 
zero  the  terms  on  the  righthand  side  of  Eq.  (A-86)  that  are  multiplied  by  p.  since  in 
the  region  of  softening  <  1  and  the  terms  P^  and  k^1  are  considerably  less  than 
unity.  Then  ^  is  given  by  Eq.  (A-81)  while  0  is  given  by  Eq.  (A-87).  i.e.. 


1  -  fl  =  jl  -  (1  -  b) 


[3  I1' 

c)2  C1 i 


(A-87) 


in  view  of  Eq.  (A-74).  Thus  a ^  may  now  be  found  upon  use  of  Eqs.  (A-80).  (A-81) 
and  (A-68)  and  is  given  below: 


A-22 


(A-88) 


2a 


»1  -  E(1  -  (1  -  *,(l  -  (1  -  b)cj|  ej' 

In  the  very  specific  case  where  m  =  b  =  1 


*1  *  E£l* 


-2(k  *  iclik 


(A-89) 


In  other  words,  directional  and  isotropic  degradation  compound  the  effect  of  softening 
in  simple  tension.  As  we  shall  see  this  is  not  the  case  in  compression.  For 
comparison  when  isotropic  damage  is  absent,  see  Eq.  (A-27b). 


(B)  SIMPLE  COMPRESSION 

In  this  case  Eqs.  (A-79)  and  (A-80)  apply  but  now  f  =  1  Thus  in  this 
specific  case 


*2f2  *  '  »  el 

Oj  =  E£j.2 


(A-90) 

(A-91) 


As  opposed  to  tension,  softening  in  compression  is  due  entirely  to  volumetric  damage. 
The  appropriate  rate  equations  in  the  presence  of  monotonic  loading  are: 


k(0)d6 


dfl  =  (1  -  0) 


■4 1  • 


i  _  k  cn~) 

V  V 


2  2) 


de. 


(A-92) 


(A-93) 


where  Eqs.  (A-71).  (A-73),  (A-90).  and  (A-92)  were  used.  Eqs.  (A-92)  and  (A-93) 
are  coupled  and  cannot  be  integrated  directly.  We  obtain,  therefore,  approximate 
solutions  of  these  equations  in  order  to  study  their  qualitative  trends.  To  that  end 
we  let  k(Q)  be  a  constant  k  and  define  c*  by  the  relation 


A-23 


(A-94) 


Equations  (A-92)  and  (A-93)  may  now  be  integrated  to  give: 


1 

*2=  {l  -  (1  -  mjkej1"1"  (A-95) 

1 

1  -  D  =  {l  -  (1  -  bjc'ej1'6  (A-96) 

1 

-»e1  =  e2f2  =  e2{l  -  (1  -  m)k  Cg}1"1"  (A-97) 

2a 

ox  =  Ee^l  -  (1  -  bjc’ej1^  (A-98) 

In  Figure  A. 5  we  show  an  experimental  relation  between  the  axial  and  transverse 
strains  in  the  case  of  a  concrete  block  under  axial  compression,  according  to  Van 
Mier  (1984).  and  the  analytical  relation  according  to  Eq.  (A-97)  when  v  =  0.35. 
k  =  70  and  m  =  2.4.  Also  shown  is  a  compressive  stress-strain  relation  for  the 


same  material  and  specimen.  The  approximate  analytical  axial  stress-axial  strain 
relation  obtained  on  the  basis^  of  Eqs.  JA-97)  and  (A-98)  is  shown  in  Figure  A. 6 
when  2a/l-b  =  1.5  and  (l-b)c  =  40  (c  assumed  constant).  For  the  same  model,  a 
comparison  between  the  predicted  and  measured  volume  strain  as  a  function  of 
the  axial  strain  is  shown  in  Figure  A. 7.  From  this  figure,  it  is  evident  that  the 
model  describes  the  observed  dilatancy  quite  well. 

A. 7.3  Discussion. 

The  extensive  and  painstakingly  precise  investigations  by  Van  Mier  into  the 
compressive  deformation  and  fracture  of  concrete  make  it  abundantly  clear  that  the 
softening  branch  of  the  axial  stress-strain  curve  represents  the  behavior  of  a  geometric 
structure  rather  than  that  of  a  uniform  material.  The  experiments  in  tension  by 
Gopalaratnam  and  Shah  (1985)  lead  to  similar  conclusions.  There  are  therefore 
significant  statistical  variations  in  the  behavior  of  concrete  in  the  "softening  regime." 

The  continuum  damage  theory  which  we  have  proposed  in  the  previous  sections 
captures,  we  believe,  many  of  the  essential  features  of  the  constitutive  behavior  of 
fracturing  concrete.  Further  aspects  of  this  approach  shall  be  pursued  in  our  future 
research. 
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Figure  A. 5 


Relation  between  axial  and  transverse  strains  under  conditions  of 
simple  compression. 
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This  part  describes  briefly  the  experimental  procedure,  the 
equipment  used,  and  results  of  the  tensile  tests  performed 
on  slotted  gray  cast  iron  specimens.  It  also  describes  the 
same  parameters  for  compression  tests  on  concrete  specimens 
(slotted  and  unslotted).  All  of  this  experimental  work  is  a 
part  of  the  ongoing  research  effort  in  the  field  of  damage 
modeling  of  brittle  materials. 

Materials : 


Two  different  materials  have  been  used  for  this  study  so  far. 
Gray  cast  iron  of  ASTM  Class  35  and  low  to  medium  strength 
concrete  were  chosen  for  this  project  based  on  their  mechanical 
properties.  Results  of  tensile  testing  of  gray  cast  iron  and 
compression  testing  of  concrete  will  be  described  in  separate 
sections  on  the  following  pages. 

TENSILE  TESTING  OF  GRAY  CAST  IRON: 


Specimen  Geometry: 


As  described  in  the  initial  report  under  this  research  contract 
(July,  1988),  the  "standard"  specimen  geometry  had  an  hour-glass 
shape  with  a  reduced  cross-section  of  0.75x0.13  in.  in  the  gauge 
section  of  the  specimen.  Figure  1  shows  schematically  the  specimen 
geometry.  Slots  of  various  lengths  and  angles  to  the  horizontal 
axis  were  cut  into  these  specimens  using  the  EDM  (Electrical 
Discharge  Machining)  process.  The  width  of  these  slots  was  kept 
constant  at  0.020  in.,  the  minimum  possible,  given  the  available 
technology.  Multiple  slots  were  also  cut  into  some  specimens  to 
determine  their  effect  on  fracture  stress  of  this  material.  These 
slotted  specimens  are  shown  schematically  in  Figs.  2  and  3. 

Figs.  4  and  5  show  pictures  of  typical  specimens  in  the  two 
categories . 


Standard  Specimen  Geometry 


Spec: 
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Fig.  3:  Slotted  Specimens  -  Angled  Slots 
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Experimental  Procedure: 


Tensile  tests  were  performed  on  all  the  specimens  using  similar 
test  parameters.  Loading  and  unloading  sequences  were  employed 
on  some  specimens  while  simple  monotonic  tensile  loading  was 
used  on  others.  Special,  large  gauge-length  strain  gauges  were 
applied  to  slotted  specimens  to  obtain  over-all  (average)  strains 
rather  than  localized  strains.  These  were  "quality-control" 
measurements  and  were  not  used  in  analytical  calculations. 

Results : 


Table  1  below  shows  the  experimental  parameters  and  results  of 
testing  performed  on  slotted  specimens.  In  all  cases,  slots  were 
horizontal  i.e;  across  the  gauge  section  of  the  specimen. 


Table  1:  Test  Results  -  Slotted  Specimens 
(Horizontal  Slots) 


Specimen  # 


Slot  Size 
(&  No.  of  Slots) 


8 

9 

10 
11 
12 


1/8  in. 
3/16  in. 
1/4  in. 
3/8  in. 
1/2  in. 


Fracture  Strength 
(psi ) 


20,833 

17,708 

17,708 

12,500 

8,542 


13 

15 


3/8  in.  (3  slots) 
3/8  in.  (5  slots) 


13,542 

12,292 
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Table  2  shows  specimen  parameters  and  fracture  strength  of 
slotted  specimens  where  slots  were  machined  at  an  angle  to  the 
horizontal  axis  as  shown  in  Fig.  3.  Slot  size  in  all  cases 
was  kept  constant  at  3/8  in. 


Table  2:  Test  Results  -  Slotted  Specimens 
(Angled  Slots) 


Specimen  #  Angle  of  Slot  Fracture  Strength 

(degrees)  (psi) 


18 

0 

13,333 

19 

0 

13,542 

20 

15 

11,094 

* 

21 

15 

17,448 

★ 

22 

30 

13,938 

23 

30 

13,583 

24 

45 

14,896 

25 

45 

14,938 

26 

60 

20,833 

* 

27 

60 

23,032 

* 

28 

75 

20,781 

29 

75 

21,553 

30 

90 

26,167 

31 

90 

28,156 

35 

60 

17,406 

@ 

Experimental 

scatter 

Extra  specimen  to 

"double-check"  the 

fracture  strength 

of  specimens 

with 

a  60  degree  slot 

Discussion : 


As  can  be  seen  from  Table  1,  the  fracture  strength  goes  down 
quite  consistently  as  the  slot  size  increases.  This  result  is 
very  much  as  expected  and  there  are  no  surprises  here.  What 
is  more  interesting  to  note  from  Table  1  is  the  insensitivity 
of  the  material  to  the  extent  of  "damage"  when  one  compares  the 
fracture  strengths  of  specimens  containing  one,  three,  or  five 
slots  of  the  same  size.  Note  that  the  specimen  failed  at  the 
same  location  (top  slot)  in  the  case  of  multiple  slots.  This 
is  an  aspect  of  fracture  and  damage  modeling  which  should  be 
explored  further  in  future  work. 
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Table  2  shows  the  fracture  strength  increasing  as  the  angle  of 
the  slot  with  the  horizontal  axis  increases.  There  is,  however, 
considerable  scatter  in  the  data  for  the  15  degree  and  60  degree 
specimens.  One  likely  reason  may  be  the  material  variability 
although,  so  far,  very  consistent  material  behavior  has  been 
observed  in  the  course  of  this  study.  This  material  variability 
may  have  arisen  from  different  temperature  histories  experienced 
by  the  different  parts  of  the  original  plate  material  during 
processing  and  heat  treatment.  So  a  certain  amount  of  scatter 
in  the  data  was  to  be  expected. 

There  is  also  a  change  in  the  fracture  mode  as  one  goes  from 
the  zero  to  45  degree  slot  group  to  45  to  90  degree  slot  group. 
In  the  first  case,  mode  I  fracture  is  the  predominant  one  while 
it  is  a  "mixed-mode"  or  mode  II  fracture  for  the  second  group. 
However,  this  change  in  fracture  mode  is  continuous  as  the  slot 
angle  changes  and  should  not  cause  variations  in  the  fracture 
strength  at  a  constant  angle.  Again,  this  is  an  extremely 
interesting  aspect  of  fracture  which  should  be  pursued  further. 
It  becomes  much  more  important  for  anisotropic  and  orthotropic 
materials  like  the  advanced  composites  where  fracture  is  almost 
always  of  a  mixed-mode  nature. 


COMPRESSION  TESTING  OF  CONCRETE: 


Various  concrete  mixes  were  tried  to  see  if  an  acceptable  mix 
could  be  arrived  at  for  production  in  a  small-batch  environment 
and  would  also  give  adequate  mechanical  properties.  The  final 
choice  of  a  typical  concrete  mix  was  as  follows: 

Cement:  3.0  lbs. 

Sand  (Fine):  7.26  lbs. 

Sand  (Coarse):  13.0  lbs. 

Water:  1.23  -  1.45  lbs. 

Water/Cement  Ratio:  0.41  -  0.48 

NOTE:  A  wet  coarse  sand  was  used  in  all  cases.  It  was  soaked 

for  15  minutes  before  use. 

Specimen  Geometry: 


Straight,  prism  type  concrete  specimens  were  used  for  this  study. 
Both  slotted  and  unslotted  specimens  were  tested  to  see  the 
effect  of  "damage"  on  fracture  strength  of  concrete.  Slots  were 
"manufactured"  into  these  specimens  by  an  ingenious  design  of 
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the  metal  mold  and  pre-placing  lubricated  metal  plates  into  the 
concrete  mix  at  the  time  of  pouring.  Molds  were  turned  over  and 
these  metal  plates  moved  at  regular  intervals  to  make  sure  that 
the  gravel  does  not  settle  due  to  gravity  at  the  bottom  of  the 
mold  and  the  plates  do  not  get  stuck  in  the  specimen  during  the 
curing  process.  The  specimen  dimensions  were  kept  constant  for 
all  tests  whether  these  were  slotted  or  unslotted  specimens. 

The  specimen  dimensions  are  given  below: 

Cross-section:  3x3  in.  square 

Height:  5  in. 

Vertical  slots  of  one,  two,  or  three  inch  lengths  were  used  in 
the  case  of  slotted  specimens.  The  slot  was  parallel  to  the  long 
axis  of  the  specimen  as  shown  in  Fig.  6. 

Experimental  Procedure: 


Compression  tests  were  performed  on  all  the  concrete  specimens 
using  similar  parameters.  A  loading  rate  of  2,000  lbs.  force 
per  minute  was  employed  in  all  cases.  recision  mechanical  dial 
gauges  were  used  to  measure  the  overall  contraction  (axial  strain) 
and  lateral  expansion  (transverse  strain)  of  the  specimens. 
Specimen  ends  were  ground  in  some  cases  to  give  flat  faces  which 
would  also  be  parallel  to  each  other  to  ensure  concentric  loading. 

Results : 


Table  3  on  the  next  page  shows  the  mechanical  properties  of 
typical  slotted  and  unslotted  specimens. 
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Table  3:  Mechanical  Properties  of  Concrete 
Specimens 


Specimen  # 

Slotted? 

Slot  Size 

Fracture  Strength 

( in.  } 

(psi ) 

H-l 

No 

1,600 

H-2 

No 

1,800 

1-1 

Yes 

1  in . 

1,541 

1-2 

Yes 

2  in . 

1,618 

1-3 

Yes 

2  in . 

1,489 

1-4 

Yes 

3  in. 

1,718 

H-3 

Yes 

3  in. 

1,756 

NOTE: 

1.  "Fracture 

Strength"  here 

means  the  maximum 

stress  level 

which  the  specimen  attained  before  the  onset  of  massive 
damage  and  crumbling. 

2.  Specimen  batches  H  and  I  were  identical  in  terms  of  their 
fabrication  specifications  and  mechanical  properties. 
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Discussion: 


It  is  quite  interesting  to  note  from  Table  3  that,  in  spite 
of  the  fact  that  specimen  batches  H  and  I  were  two  different 
groups  of  specimens,  the  fracture  strength  in  all  cases  stays 
relatively  constant.  In  other  words,  the  material  does  not 
see  the  slots,  whether  small  or  large,  as  "damage"  and  behaves 
as  if  these  slots  were  not  even  there!  This  point  is  discussed 
at  length  in  the  main  body  of  the  report. 

The  micro-cracks  which  developed  in  all  specimens  had  the  same 
pattern  of  development  as  the  load  increased.  Most  predominant 
orientation  of  these  cracks  was  along  the  loading  axis  i.e; 
cracking  due  to  lateral  expansion  of  the  specimen.  This  mode 
of  vertical  splitting  was  seen  in  all  cases  and  the  presence 
of  vertical  slots  did  not  seem  to  affect  it  in  any  major  way. 
The  damage  did  not,  for  example,  start  at  the  ends  of  these 
slots.  Fig.  7  shows  this  vertical  splitting  in  some  slotted 
specimens. 

All  of  the  above  observations  are  extremely  interesting  from 
the  point  of  view  of  the  definition  of  "damage"  in  a  given 
material.  When,  for  example,  a  crack,  notch,  or  a  geometrical 
discontinuity  becomes  "damage"  and  is  critical  for  the  service 
life  of  a  material/structure?  These  questions  should  be  looked 
into  very  carefully  through  future  work  in  the  area  of  fracture 
control  and  damage  tolerance. 


tu  t-3  CU 


Slotted  Concrete  Specimens  showing 


